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ABSTRACT 


Translocation  of  the  Global  Positioning  System  has  proved  to  be  a  highly 
accurate  method  of  position  determination  for  onshore  and  airborne  navi¬ 
gation,  but  it  had  not  been  previously  evaluated  for  nearshore  hydro- 
graohic  surveys.  The  technique  of  translocation  for  hydrographic  opera¬ 
tions  involves  the  simultaneous  reception  of  signals  from  the  GPS 
satellites  by  two  independent  receivers;  one  receiver  onboard  the  survey 
vessel  and  one  located  at  a  known  stationary  position.  A  position 
correction  is  obtained  at  the  stationary  receiver,  which  is  then  applied 
to  the  shipboard  receiver  online  or  during  postprocessing.  An  accuracy 
determination  of  translocated  GPS  was  conducted  at  the  Naval  Postgraduate 
School,  Monterey,  CA,  May,  1981.  Two  methods  of  positioning  were  used 
for  comparison  with  GPS:  1)  a  least  squares  solution  of  three  lines  of 
position  observed  from  three  Wild  T-2  theodolites,  and  2)  a  position 
determined  from  ARTEMIS,  a  range-azimuth  short  range  microwave  position¬ 
ing  system.  Translocated  GPS  accuracies  of  10  meters  were  determined. 

It  is  anticipated  that  greater  accuracies  will  be  obtained  by  using  a 
more  sophisticated  receiver  and  more  advanced  processing  methods. 
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Shipboard  operations  for  the  purpose  of  oceanographic 
research,  coaaercial  transport,  tactical  operations,  or 
hydrographic  surveying  require  accurate  and  dependable 
positioning  systeas.  The  operation  of  aost  electronic 
positioning  systeas  is  haapered  by  ataospheric  effects  on 
the  propagated  signal  and  design  liaitations  which  restrict 
operating  ranges,  positional  accuracy  and  dependability.  & 
universal  positioning  systea  that  provides  accurate 
inforaation  to  the  user  regardless  of  ataospheric  conditions 
and  that  is  free  of  range  liaitations  is  desirable.  with 
the  advent  of  satellite  positioning  techniques,  such  a 
systea  is  now  available. 

The  Navy  Navigation  Satellite  Systea  (NAVSAT) ,  or 
TRANSIT,  has  been  operational  since  1964  and  has  proved 
dependable  and  sufficiently  accurate  for  certain  users. 
Although  developed  for  ailitary  uses,  the  application  of 
TRANSIT  in  the  civilian  coaaunity  has  continued  to  expand  at 
an  exponential  rate.  The  systea,  however,  is  liaited  by  the 
tiae  interval  between  fixes,  which  is  90  ainutes  on  the 
average,  and  the  requireaent  for  precise  knowledge  of  the 
ship's  velocity.  Although  the  systea  does  not  provide 


One  aethod  of  improving  system  accuracy  is  by  using  two 
receivers  simultaneously  in  a  technique  termed 
translocation,  or  differential  mode.  The  method  is  based  on 
the  assumption  that  signal  propagation  errors  to  two 
proximate  receivers  are  nearly  the  same.  One  receiver  is 
placed  over  a  known  position,  and  the  variation  of  the 
observed  position  with  time  is  determined  relative  to  the 
absolute  position.  This  variation  is  applied  to  the 
observed  positions  determined  by  a  simultaneously  operating 
receiver  located  at  another  location — on  a  survey  vessel, 
for  example.  The  assumption  underlying  this  technique  is 
that  the  distance  between  the  two  units  is  not  sufficient  to 
introduce  a  significant  additional  propagation  error.  The 
resulting  position  of  the  translocated  receiver  should  be 
more  accurate  than  the  position  obtained  using  a  single 
receiver. 

For  hydrographic  operations,  the  translocated  Global 
Positioning  System  (GPS)  provides  an  alternative  to 
presently  used  positioning  systems.  The  satellites  orbit  at 
20,000  km,  which  theoretically  would  allow  receiver 
separation  of  as  much  as  500  km  without  causing  significant 
additional  system  error.  frith  this  flexibility,  one 
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sufficient  positional  accuracy  for  hydrographic  surveying, 
* 

it  has  been  used  for  geodetic  application  by  the  Defense 
Napping  Agency  (DMA)  using  fixed  shore-based  receivers  to 
establish  positioning  networks  in  inaccessible  areas 
worldwide.  [Ref.  1] 

The  Global  Positioning  Systea,  or  NAVSTAR,  is  a  more 
sophisticated  systea  designed  to  replace  TRANSIT  in  the 
aid-1 980s.  Six  satellites  are  currently  in  orbit  with 
twelve  aore  to  be  launched  before  the  operational  date  of 
1987.  The  18-satellite  orbit  configuration  will  provide  a 
aulti-user,  passive  systea  with  24-hour  availability,  global 
coverage,  and  possible  fix  intervals  of  less  than  one 
second.  The  operational  network  consists  of  three 
eleaents — the  satellites,  the  ground-based  control  segaent, 
operated  by  the  Space  and  Nissile  Systens  Organization 
(SANSO) ,  and  the  user*s  receiver.  Prototype  receivers  have 
been  tested  under  optiaua  conditions  by  SANSO  at  the  Tuna 
Proving  Ground,  Arizona,  on  both  aobile  and  stationary 
platforas.  Results  show  excellent  systea  stability  and 
accuracy  when  conpared  to  currently  used  positioning 
systeas.  [Ref.  2] 
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shorebased  receiver  could  be  placed  at  a  central  location 
and  a  survey  of  several  hundred  kilometers  of  coastline  or 
offshore  area  could  be  coapleted  with  only  one  additional 
receiver  on  the  survey  vessel.  [Ref.  3]  The  value  of  this 
flexibility  would  be  realized  in  the  significant  decrease  in 
survey  planning,  cost,  and  time  compared  to  that  spent 
presently  in  conducting  hydrographic  operations. 

The  accuracy  of  the  systea  would  be  sufficient  for  all 
survey  scales.  The  Rational  Ocean  Survey  (MOS)  defines  the 
required  accuracy  of  a  positioning  systea  to  be  no  more  than 
one-third  of  the  survey  accuracy,  which  is  1.5  ailliaeters 
at  the  scale  of  the  survey.  [Ref.  4]  Mon- translocated  GPS 
has  been  shown  to  be  of  sufficient  accuracy  for  1:80,000 
scale  surveys  and  saaller.  [Ref.  5] 

The  object  of  our  research  is  to  determine  whether  GPS 
at  the  present  operational  level  is  adequate  for  1:10,000 
scale  surveys.  To  do  this,  a  position  determined  by  two  GPS 
receivers  in  the  translocated  aode  is  coapared  to  the  saae 
position  determined  by  two  other  aethods.  If  translocated 
GPS  can  be  shown  to  achieve  specified  inshore  survey 
accuracies,  the  systea  would  provide  sufficient  accuracies 
at  all  saaller  scale  surveys. 
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Part  of  th«  difficulty  in  determining  this  level  of 
accuracy  concerns  the  reference  system  used.  (lost  geodetic 
positions  published  by  the  National  Geodetic  Survey  (NGS) 
are  based  on  the  North  American  Reference  Datua  of  1927 
(NAD-27)  and  the  Clarice  1866  Ellipsoid,  the  best 
approximation  to  the  shape  of  the  ear*'1,  f ■  the  continental 
United  States.  The  satellite  s yst**** ,,  k-*«ever,  use  a 
reference  systea  which  is  not  bast:'1,  mv  portion  of  the 
earth's  surface.  A  mass-centered  ei«kys©id  referenced  to 
the  earth's  center  of  mass  has  be.'}*,  computed  using  observed 
gravity  and  astronomical  data.  The  best  fit  to  this  global 
ellipsoid  is  the  World  Geodetic  Systea  of  1972  (WGS-72) , 
which  is  utilized  by  both  *  TRANSIT  and  GPS.  Transformation 
equations  exist  which  convert  NAD-27  coordinates  to  wgs-72 
coordinates,  but  uncertainty  in  the  gravity  model  produces 
errors  as  large  as  10  aeters  due  to  the  transformation 
process  alone.  [Ref.  6] 

To  establish  all  positions  on  the  saae  reference 
ellipsoid  for  this  research,  a  TRANSIT  receiver  was  acquired 
from  DNA  and  a  first  order  Doppler  station  was  established 
at  the  central  location  of  the  test  area.  A  third  order 
geodetic  survey  was  conducted  froa  the  TRANSIT  station  to 
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establish  two  additional  stations  on  the  BGS-72  ellipsoid. 
Two  single  channel  Manpack  GPS  receivers  were  obtained  fros 
SAfiso.  Stability  testing  and  a  receiver  cosparison  test 
were  conducted  prior  to  shipboard  testing  to  deteraine  a 
quantitative  aeasure  of  receiver  perfocaance.  One  of  the 
receivers  was  then  placed  over  the  TRANSIT  station  aark  and 
the  other  aboard  a  survey  vessel  which  naneuvered  along 
predeterained  track  lines.  Data  were  recorded  via  an 
interface  to  an  HP-9825  coapater  and  printer  at  each  site. 
Siaultaneously,  two  additional  independent  aeasureaents  were 
aade  of  ship  position — one  using  a  range-aziauth  short 
range  aicrowave  positioning  systea  (ARTEMIS)  and  the  other  a 
position  deterained  by  a  least  squares  solution  of  three 
lines  of  position  observed  froa  three  1-second  theodolites 
set  on  the  three  pre-established  BGS-72  locations.  The 
coaparison  of  the  ARTEMIS  and  theodolite-derived  ship 
positions  with  the  translocated  GPS  ship  position  provided 
the  test  results. 
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II.  BQOIPMENT  DESCRIPTION  A ND  OPERATIONAL  CAPABILITIES 

Two  positioning  systems  were  chosen  for  comparison  with 
the  Global  Positioning  System:  a  theodolite  network  and  a 
microwave  positioning  system.  The  theodolite  network 
consisted  of  three  shore-based  units  so  as  to  provide  three 
lines  of  position.  The  microwave  positioning  systen 
provided  a  range  and  an  azimuth  from  a  known  location  on 
shore.  Both  systems  are  accurate  at  short  ranges  and 
represent  two  contrasting  methods  of  positioning — visual  and 
electronic.  The  theodolite  method  of  positioning  at  short 
range  has  been  and  continues  to  be  used  for  hydrographic 
operations  as  an  alternative  to  electronic  systems  when 
gecmetry  and  system  accuracy  limitations  impose  the  need  for 
a  more  flexible  positioning  system.  The  requirement  for 
experienced  theodolite  observers  as  well  as  the 
communication  and  logistics  difficulties  involved  limit  the 
use  of  this  method  to  special  cases. 

This  chapter  describes  the  AETENIS  short-range 
electronic  positioning  system,  the  operation  of  the 
theodolite  network  as  used  during  the  comparison  test,  and 
the  TRANSIT  satellite  system. 
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A.  ARTEMIS  POSITIONING  SYSTEM 


The  ARTEMIS  position  fixing  system  is  a  short-range 
microwave  systea  (to  30  ka) ,  built  by  Christiaan 
Huygenslaboratoriua  B.  V . ,  in  the  Netherlands.  It  has  been 
used  operationally  (priaarily  in  Europe)  since  1972,  when 
the  systea  was  introduced  at  the  International  Hydrographic 
Conference  in  Monaco.  At  the  present  tine,  fifty  systeas 
are  in  use.  The  systea  is  unique  in  design,  is  highly 
accurate  at  short  ranges,  and  is  particularly  useful  for 
snail  harbor  surveys,  river  surveys  and  relative  positioning 
uses.  Por  this  study,  the  ARTEMIS  systea,  along  with 
technical  assistance,  was  obtained  froa  the  American 
distributor,  MAHINAV  Corporation  of  Houston,  TEXAS. 

ARTEMIS  is  a  range-azimuth  positioning  system  which 
eaploys  a  unique  tracking  method.  A  stationary  or  fixed 
unit  is  placed  over  a  known  position  on  shore.  Using  a 
telescope  that  is  mounted  on  the  antenna,  the  antenna  is 
sighted  onto  another  known  point  and  the  known  azimuth 
between  the  two  locations  is  entered  into  the  systea  via  a 
digital  display  on  the  fixed  unit.  The  unit  is  then  locked 
into  place  for  operation.  A  mobile  unit  is  placed  on  the 
vessel  whose  position  is  desired. 
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The  mobile  and  fixed  antennas  resemble  radar  antennas 
which,  while  operating,  track  each  other  so  that  the  two 
antennas  are  always  parallel  to  each  other  and  perpendicular 
to  the  line  of  position  between  them.  A  microwave  link  is 
established  between  the  two  antennas  during  periods  of 
measurement  and  information  is  relayed  between  the  mobile 
and  fixed  units.  (A  voice  link  is  also  available  during 
operation.)  The  fixed  unit  provides  the  direction  and  the 
azimuth  of  the  mobile  unit  with  respect  to  the  point  sighted 
on  by  the  telescope.  At  the  same  time,  the  mobile  unit 
transnits  a  coded  interruption  of  the  signal  to  the  fixed 
unit,  which  replies  with  the  same  interruption.  The  time 
difference  between  the  transmission  and  reception  at  the 
aobile  unit  determines  the  distance  between  the  two 
antennas. 

Distance  accuracy  is  determined  by  the  number  of 
measurements  averaged  over  tine  and  is  available  in  two 
modes — a  dynamic  mode  which  averages  1000  measurements  or  a 
static  mode  which  measures  10,000  measurements  per  displayed 
distance.  The  azimuth  and  distance  information  is  displayed 
by  the  mobile  unit  so  that  personnel  are  not  required  at  the 
stationary  unit  once  operation  has  begun.  The  system  can  be 
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operated  in  long  or  short-range  nodes,  which  deteraine  the 
aaount  of  power  output  for  the  system. 

The  microwave  frequencies  used  are  9.2  and  9.3  GHz. 

Problems  arising  with  signal  interference,  observed  with  all 
microwave  positioning  systems,  therefore  can  occur  with 
ARTSMIS.  The  regions  of  range  holes,  or  multipath 

interference,  can  be  computed  for  up  to  10  km  range  by  prior 
knowledge  of  the  antenna  elevations  above  sea  level  and  the 
operating  distance  from  the  shore- based  antenna.  [Ref.  7] 

Each  antenna  is  designed  with  a  22°  vertical  beam  width 
and  a  2°  horizontal  beam  width.  A  power  source  of  22  volts 
DC  per  unit  is  required,  with  an  average  current  load  of  2.5 
amps. 

The  accuracy  of  the  system  depends  on  both  the  distance 
and  the  angle  measured.  The  angle  measurement  has  a  2  a 
probability  of  measuring  within  2’  of  arc  (0.033°  >  .  The 

distance  measurement  has  a  2  G  probability  that  .he  mean 

distance  value  will  vary  within  +  1.5  meters. 

One  of  the  advantages  of  ARTEHIS  over  other  electronic 
systems  is  that  the  two  lines  of  position  defining  the  fix 
always  intersect  at  90°  .  A  unique  advantage  is  that  the 
coverage  around  the  fixed  unit  can  be  circular,  with  a 


aaxiaua  range  of  operation  to  30  ka,  dependent  on  antenna 
heights.  The  short  range  applications  of  the  ABTEMIS  fill  a 
need  not  provided  by  aany  currently  used  positioning 
systeas. 

Por  this  study,  the  ABTESIS  was  used  in  conjunction  with 
the  positions  derived  fron  the  theodolite  network  to 
deteraine  the  best  estiaate  of  the  ship*s  position.  The 
accuracy  of  the  ABTEMIS  was  more  than  adequate  for  the 
ranges  encountered  during  the  test,  and  the  ease  of  use  and 
operational  versatility  of  the  systea  made  its  use  very 
advantageous. 

B.  THEODOLITE  HETWORK 

A  network  of  three  Wild  T-2  theodolites  was  used  as  the 
control  positioning  systea.  The  theodolite  is  widely  used 
for  surveying  as  a  leans  of  aeasuring  angles  between  two 
observational  locations — either  in  a  horizontal  plane  for 
traverse  or  triangulation  surveys,  or  a  vertical  plane  for 
astronoaical  observations. 

The  use  of  three  theodolites  as  a  short-range 
positioning  systea  was  an  accepted  aethod  of  positioning 
prior  to  developaent  of  electronic  distance  aeasuring 
systeas.  There  reaain  disadvantages  when  using  electronic 
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positioning  systems  for  hydrographic  operations  at  short 
ranges  (1-2  km).  Jlost  short-range  electronic  systeas  are 
designed  for  measurement  of  position  at  line-of-sight 
ranges.  This  region  includes  a  large  portion  of  inshore 
hydrographic  surveys.  Inshore  positioning  involves  aore 
thorough  planning  to  obtain  the  required  accuracy  for  the 
scale  of  the  survey.  Por  example,  in  a  harbor  or  cove  it 
may  be  difficult  to  place  two  range-range  electronic  units 
so  that  the  area  to  be  surveyed  is  not  within  the  region  of 
positioning  systea  uncertainty,  or  30°  (radially)  on  each 
side  of  the  baseline  between  the  two  units.  To  reposition 
the  electronic  units  as  often  as  is  required  to  obtain 
optiaua  geoaetry  is  tiae-consuaing  and  requires  manpower  and 
extensive  planning,  and  can  result  in  operational  delays. 
Interference  of  the  signal  due  to  reflection  from  a  calm  sea 
surface,  or  aultipathing,  is  also  a  problem  at  short  ranges 
for  aicrowave  systeas. 

Theodolites  can  be  set  at  any  intervisible  location  with 
preplanned  geometry  and  can  be  used  for  the  entire  region  to 
be  surveyed.  This  aethod  requires  less  expensive  equipment 
(three  theodolites) ,  more  manpower  (a  recorder-  and  an 
observer  at  each  site) ,  communication  between  the  three 


locations  (via  portable  radio,  for  example)  and. reasonable 
veather  conditions.  The  basis  of  this  procedure  is  that 
three  angles  observed  from  known  locations  to  a 
predeterained  point  on  the  vessel  constitute  three  lines  of 
position.  After  removing  all  possible  systeaatic  errors  and 
through  a  least  squares  solution  of  the  three  lines  of 
position,  the  position  of  the  vessel  is  deterained.  A  major 
disadvantage  is  that  observational  errors  are  not  deterained 
in  real  time  because  the  data  cannot  be  processed 
immediately. 

The  basic  procedure  for  the  method  requires 
predeterained  sites  which  are  known  relative  to  each  other. 
This  may  be  achieved  by  using  published  geographic  positions 
or  by  conducting  a  local  plane  survey  to  establish  new 
positions  using  preestablished  stations  nearby.  Either  way, 
all  three  locations  must  be  based  on  the  same  datum  and 
preferably  permanently  positioned  by  a  disk  or  a  station 
mark. 

Each  theodolite  was  aligned  using  a  plumb  bob  or  optical 
plummet  over  one  of  the  known  station  marks  or  disks.  By  a 
predetermined  scheme,  each  site  was  marked  so  as  to  be 
clearly  visible  to  the  other  sites  (lights  were  used  for 
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this  study)  .  The  telescope  of  each  theodolite  was  sighted 
onto  one  of  the  other  sites  and  the  horizontal  circle  set  to 
read  an  arbitrary  value  close  to  0°.  The  crosshairs  in  the 
telescope  were  set  to  intersect  the  mark.  When  the 
theodolite  was  exactly  sighted,  the  horizontal  circle  was 
read  and  recorded.  (At  this  tine  the  lights  were  turned 
off.)  The  operation  began  when  all  three  sites  coapleted 
this  step,  called  the  "initialling  procedure",  or  "initial 
pointing".  The  instruaent  was  sighted  on  a  point  on  the 
survey  vessel  that  was  visible  to  all  operators.  The  ship 
was  tracked  by  the  operators  by  using  the  aicroaeter  wheels 
on  the  instruaent  until  a  "aark"  was  called  over  the  radio 
by  the  recorder  at  the  control  theodolite  location.  At  that 
instant,  the  operator  stopped  tracking  and  the  horizontal 
circle  was  read  and  recorded.  This  observation  constituted 
a  fix.  After  a  series  of  fixes  were  recorded,  or  when  the 
survey  vessel  reached  a  predeterained  location,  a  set  of 
fixes  was  coapleted.  The  lights  were  then  reactivated  at 
each  site,  each  observer  intersected  the  original  aark,  the 
horizontal  circle  was  again  read  and  the  "final  pointing" 
was  recorded.  If  necessary,  the  instruaent  was  relevelled. 
The  process  was  then  repeated  by  starting  with  the 
initialling  procedure. 
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It  was  assumed  that  the  theodolites  were  in  good 
operating  condition.  Collimation  error  had  previously  been 
determined.  The  observers  were  kept  at  the  same  locations 
throughout  the  test  to  minimize  the  error  due  to  individual 
observer  variations.  Since  only  one  portion  of  the 
horizontal  circle  was  used,  and  no  reverse  readings 
observed,  a  systematic  error  was  introduced.  This  was 
assumed  to  be  minimized  during  computation  of  the  vessel 
position. 

The  observations  were  corrected  for  collimation  error 
and  for  pointing  error  by  evenly  distributing  the  difference 
between  the  initial  and  final  pointing  values.  The 
resulting  angles  were  processed  using  an  intersection 
computation  with  a  least  squares  fit  applied  to  optimize  the 
result.  The  computations  for  this  study  were  processed 
using  an  HP-9815  computer  and  software  prepared  by  Cdr 
Ludvig  Pfeifer,  NOAA.  The  resulting  position  was  used  as 
the  best  estimate  of  the  position  of  the  ship. 

C ,  TRANSIT  POSITIONING  SYST2H 

The  TRANSIT,  or  Navy  Navigation  Satellite  System,  is  the 
most  widely  utilized  satellite  positioning  system  at  the 
present  time.  The  network  of  satellites  and  ground  tracking 


stations  has  been  in  operation  since  1964  and  is  becoming 
widely  used  by  both  the  military  and  civilian  coeaunities. 
It  is  expected  to  continue  operation  until  at  least  1995, 
when  the  Global  Positioning  Systee  is  expected  to  replace 
it. 

The  systea  consists  of  a  series  of  T3ANSIT  satellites 
(fig.  2.1.)  in  fixed  polar  orbits  of  1095  lea  elevation, 
feraing  a  network  within  which  the  earth  rotates  (fig. 2. 2.). 
The  orbits  are  spaced  so  that  satellite  passes  occur,  for  a 
given  ground  location,  at  intervals  of  between  35  and  100 
ainutes.  Each  satellite  travels  an  observable  distance  of 
between  4400  and  7000  km  per  pass,  providing  a  sufficient 
tiae  interval  and  baseline  over  which  the  user  can  record 
data. 

The  idea  of  positioning  an  observer  on  the  earth's 
surface  by  using  the  satellite  as  a  reference  system 
originated  froa  observation  of  Sputnik  1,  launched  by  the 
Soviet  Onion  in  1957.  By  recording  the  change  in  frequency 
of  the  satellite  signal  as  the  satellite  passed  overhead,  a 
Doppler  shift  could  be  measured.  The  Doppler  shift  gives  a 
spatial  relationship  as  a  function  of  tiae  between  the 
observer  and  the  satellite.  rf  the  position  of  the 
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satellite  is  known  at  the  tine  of  observation/ 
observer's  position  can  be  deterained. 

In  the  TRANSIT  system,  the  satellite  positions  are 
precisely  deterained  by  tracking  stations  located  at 
Prospect  Harbor,  Maine;  Rosaaount,  Minnesota;  and  Hahiawa, 
Hawaii.  The  tracking  information  (in  the  fora  of  the 
Doppler  frequency  shift  as  a  function  of  tine)  is  relayed  to 
the  coaputing  center  at  the  Naval  Astronautics  Group 
Headquarters,  Pt.  Mugu,  California.  Here  the  actual 
satellite  orbits  are  coaputed.  This  inforaation  is  used  to 
predict  orbit  inforaation.  Parameters  of  this  predicted 
orbit  are  relayed  to  two  injection  stations  at  Roseaount  and 
Pt.  Mugu.  During  the  next  satellite  pass,  the  information  is 
relayed  to  the  satellite  as  part  of  the  navigation  message, 
with  new  updates  occurring  every  12  hours. 

The  satellites  therefore  contain  predicted  orbit 
inforaation  which  is  received  by  the  user  via  a  satellite 
receiver.  Several  receivers  are  available  with  complexity 
and  positional  accuracy  dependent  upon  the  needs  of  the 
user.  A  computer  is  noraally  required  to  deteraine  the 
user's  aobile  or  fixed  position  in  real  tiae  from  a  single 
satellite  pass. 
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Por  many  users  having  a  dual  channel  receiver  and  known 
velocity,  this  aethod  of  positioning  provides  sufficient 
accuracy  with  typical  error  of  27  to  37  aeters  root  sua 
sguare.  Optiaua  accuracy  is  obtained,  however,  if  the 
receiver  position  is  fixed,  and  several  consecutive  passes 
of  the  saae  satellite  are  recorded.  During  each  pass, 
lasting  froa  10  to  16  ainutes  as  the  satellite  travels  froa 
one  horizon  to  the  other,  aany  Doppler  counts  (observations 
of  shift  in  frequency  as  a  function  of  tiae)  are  recorded. 
After  a  predeternined  nuaber  of  passes  are  stored,  the  data 
can  be  forwarded  to  the  Defense  Mapping  Agency  in  Brookaont, 
Maryland.  There  the  precise  epheaeris  data  froa  the  ground 
tracking  systea  is  applied  to  the  data  froa  the  recorded 
passes  to  provide  the  optiaua  solution.  The  resulting 
position  has  a  typical  error  of  6.3  aeters  root  aean  square 
for  a  series  of  passes,  with  a  3-diaensional  result  of  1.5 
aeters  per  axis  repeatability  after  25  precise  epheaeris 
passes.  This  aethod  is  used  to  deteraine  positions  in 
reaote  areas  by  DMA  geodetic  survey  groups  using  the 
Magnavox  AN/PRR-14  Seoceiver  (fig. 2. 3.).  [Ref.  8] 

The  TRANSIT  systea  uses  the  World  Seodetic  Systea  of 
1972  (WGS-72)  reference  ellipsoid  which  is  based  on  the 
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Universal  space  Rectangular  coordinate  Systea — an  earth, 
aass-centered  coordinate  systea.  The  BGS-72  ellipsoid 
describes  the  best  fit  of  the  entire  earth  to  a  aatheaatical 
aodel.  The  satellite  orbits  are  coaputed  using  this 
geopotential  aodel,  the  accuracy  of  which  is  based  on 
gravitational  and  astronomical  data. 

For  the  present  study,  a  Hagnavox  AN/PRR-14  Geoceiver 
was  acquired  froa  the  Department  of  Satellite  Geophysics, 
Satellite  Tracking  Branch,  DMA.  The  receiver  was  placed 
over  a  fixed  point  and  40  passes  of  one  satellite  were 
observed.  The  procedures  described  previously  for  an 
optiaua  solution  of  point  positioning  were  followed.  The 
resulting  latitude,  longitude,  and  height  relative  to  the 
WGS-72  ellipsoid  are  coordinates  of  a  first  order  TRANSIT 
position  as  defined  by  DBA  standards.  This  position  was 
used  as  the  basis  for  horizontal  control  for  the  study  (see 
test  procedures) . 
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A.  SYSTEM  DESCRIPTION 

NAVSTAR,  or  the  Global  Positioning  System,  was  designed 
for  maximum  operational  capability — high  accuracy, 
twenty-four  hour  access,  passive  user,  and  world-wide 
coverage.  A  brief  description  of  the  GPS,  and  the 
satellites  (NAVSTARs),  will  serve  as  a  background  for  the 
more  in-depth  discussion  immediately  following. 

GPS  navigation  is  based  on  the  measurement  of  four 
ranges,  or  pseudo- ranges,  from  each  of  four  satellites  with 
known,  or  predicted,  positions.  Each  range  determination  is 
computed  using  the  propagation  velocity  of  the  signal 
multiplied  by  the  satellite-to-user  travel  time.  Satellite 
transmission  time  is  included  as  part  of  the  navigation 
signal  message.  Prom  the  four  known  ranges,  the  user 
position  and  user  clock  error  are  computed  on  a  three-axis, 
orthogonal  coordinate  system.  This  computation  is  then  used 
to  transform  the  user  position  onto  the  'dGS-72  ellipsoid. 

The  remainder  of  this  chapter  will  provide  the 
background  for  GPS  development  and  the  three  major  system 
segments  of  the  GPS:  the  space  system  segment,  the  control 


segment,  and  the  user  system  segment.  Included  in  the  user 
system  segment  will  be  a  discussion  of  the  Nanpack  receiver 
as  developed  by  Texas  Instruments  Corporation  used  for  this 
study. 

With  the  successful  operation  of  the  TRANSIT  system  in 
the  early  1960's,  the  Department  of  Defense  (DOD)  began  to 
consider  the  specifications  of,  and  the  technology  required 
for,  the  development  and  implementation  of  a  second 
generation  system.  On  17  April,  1973,  the  Deputy  Secretary 
of  Defense  stated  that  a  Defense  Navigation  Satellite  System 
was  to  be  developed.  [Sef.  9]  Using  the  new  synthetic 
oscillator  technology,  and  the  Department  of  Defense  surface 
and  nearsurface  navigation  objectives  and  requirements,  a 
DOD  interdepartmental  task  force  developed  the  GPS  as  the 
second  generation  satellite  navigation  system.  With  the 
United  States  Air  Force  (USAF)  as  the  managing  agency,  the 
GPS  became  a  joint  service  program,  with  representatives 
from  the  DMA,  Army,  Navy,  Marines,  Coast  Guard,  and  the 
North  Atlantic  Treaty  Organization.  [Ref.  10]  The  Joint 
Programs  Office  was  established  within  the  Air  Force  Space 
and  Missile  Systems  Organization  to  design,  develop,  and 
implement  the  GPS  program. 


GPS  was  conceived  and  designed  to  provide  continuous. 


real-time,  world-wide  navigation  coverage  for  subsurface, 
surface  and  nearsurface  operational  vessels,  vehicles,  and 
aircraft.  Accuracy  of  better  than  10  meters  was  desired 
under  all  conditions,  with  navigational  capabilities  to 
update  and  interface  with  other  forms  of  existing  and 
projected  navigational  systems.  Another  requirement  of  the 

i 

GPS  satellites  was  nearly  trouble-free  and  simple  user 
operation,  even  in  low  signal-to-noise  environments.  An 
encoded  signal  structure  was  required  for  the  satellite 
system  for  security  reasons.  In  addition,  the  signal 
transmitted  by  the  NAVSTAB  satellites  had  to  include  other 
operational  information  in  conjunction  with  the  navigational 
message.  To  provide  for  navigational  requirements  in  event 
of  hostilities,  the  GPS  was  designed  to  permit  the  use  of 
redundant  satellites. 

To  achieve  these  objectives  and  requirements,  the  GPS 
program  was  organized  according  to  three  developmental 
tasks:  concept  validation,  full-scale  engineering 
development  and  systems  tests,  and  production  and  operation. 
Starting  in  1974,  the  concept  validation  task  was  initiated 
with  the  merger  of  the  TISATION  series  satellite  program  and 


the  US  AP  System  21 B  (a  highly  accurate  three-dimensional 
navigation  systea)  prograa.  Tests  with  the  Navigation 
Technology  Satellite,  NTS-1,  utilizing  two  rubidiua  clocks 
(accuracy  5-10  parts  per  1013)  and  an  orbit  7500  nautical 
ailes  (na)  above  the  surface  and  three  ground  based 
receivers  (called  inverted  range)  provided  concept 
verification.  [Ref,  11]  Additional  tests  were  conducted 
with  an  improved  NTS-2  satellite  positioned  at  an  altitude 
of  10,980  na  and  containing  two  cesium  clocks  (accuracies  on 
the  order  of  1  part  per  10  n).  [Ref.  12]  Orbital  height  of 
10,980  nm  permitted  the  evaluation  of  the  effect  of  the 
planned  GPS  constellation  altitude  on  navigational 
accuracies  using  the  improved  cesium  frequency  standards. 
In  1980,  NTS-3  was  launched  and  the  Advanced  Developmental 
Nodel  (ADH)  hydrogen  maser  frequency  standard  evaluated  for 
use  on  the  GPS  (providing  an  expected  frequency  accuracy  of 
1  part  per  1013).  [Ref.  13]  The  results  obtained  from  the 
NTS-3  satellite  have  not  yet  been  released.  Since  1978, 
additional  satellites  have  been  orbited  to  provide  the 
six-satellite  constellation  required  to  fully  evaluate  the 
first  developaent  task — concept  validation.  Phase  I. 
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Daring  Phase  I,  user  requirements  were  consolidated  to 
allow  the  development  and  construction  of  four  basic  types 
of  test  receivers:  continuous  tracking  receivers, 
sequential  tracking  receivers,  lew  cost  receivers,  and  the 
Hanpack  receivers.  [Ref.  14]  These  receiver  types  were 
then  constructed  for  Phase  r  testing  by  various  contractors 
and  their  perforaances  evaluated  aboard  helicopters, 
aircraft,  ships,  vehicles  and  personnel  at  the  Joint 
Prograas  Office's  testing  facilities.  Phase  I,  concept 
validation,  was  completed  in  late  1979  and  Phase  II  testing 
was  begun  in  early  1980  with  full  scale  engineering, 
developmental,  and  systems  tests. 

Phase  II  tasks,  finishing  the  operational  and  control 
segments  of  the  testing,  are  scheduled  for  completion  in 
aid-1983.  The  operational  control  station  will  be  relocated 
at  a  central  continental  united  states  site  during  this 
phase.  Additionally,  all  satellite  nonitor  receivers  at  the 
remote  aonitoring  stations  will  be  updated  to  operational 
status.  User  equip aent  and  receiver  selection  will  be  made 
in  two  increments.  Four  contractors  will  be  asked  to 
coapete  for  the  second  increaent  of  the  receiver  selection 
process  by  providing  upgraded  and  refined  user  receivers  for 


extensive  testing 

and 

evaluation. 

Two  of  the 

four 

contractors  will 

then 

be  selected 

to  continue 

the 

competition  and  will  provide  additional  receiver  design, 
refinement  and  modification,  leading  to  prototype  production 
and  extensive  field  testing. 

When  Phase  II  is  completed  in  aid-1983,  Phase  III, 
production,  will  commence  with  the  selection  of  one  of  the 
two  contractors  to  produce  the  user  receivers.  Manufacture 
of  these  receivers  will  commence  shortly  after  selection, 
and  commercial  availability  of  receivers  is  expected  in 
1985.  [Ref.  15]  The  GPS  should  be  fully  operational  with 
world-wide  coverage  using  a  constellation  containing  a  total 
of  eighteen  satellites  in  1988.  [Ref.  16,  173 

The  GPS  is  considered  to  be  composed  of  three 
inter-related  system  segments:  the  space  system  segment, 
the  control  segment,  and  the  user  system  segment.  Of  these, 
the  space  system  segment  may  be  considered  the  most 
important  since  the  other  two  reguire  the  space  segment  to 
be  operational  at  all  times.  Operationally,  the  space 
system  segment  consists  of  18  satellites  placed  in  orbits  at 
altitudes  of  10,898  nm  and  having  orbital  periods  of  about 
12  hours.  [a«f.  18] 
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The  original  plan  for  the  satellite  constellation 
consisted  of  24  satellites  in  three  planes  of  8  satellites 
each.  These  three  planes  were  to  be  placed  45°  of  longitude 
apart  and  at  an  inclination  of  63  0  to  the  equator, 
f Hef .  19]  As  a  result  of  1980  budget  cutbacks  the  space 
systen  segment  has  been  reduced  to  a  total  of  18  satellites. 
£Ref.  20]  Several  different  types  of  18  satellite 
constellations  have  been  proposed,  but  all  have  drawbacks. 
One  which  appears  the  aost  promising  is  the  ' Nonuniforn  18* 
constellation  containing  6  planes,  each  of  which  contains  3 
satellites  spaced  120°  apart  (fig. 3.1.).  The  areas  of  user 
outages  (fig. 3.2.)  are  a  disadvantage  of  this  configuration. 
The  user  outages  are  indicated  by  the  darker  areas  shown  on 
?igure  3.2.  These  outage  areas  will  appear  to  rotate  about 
the  globe  at  the  indicated  latitudes,  and  positioning  with 
the  GPS  will  be  intermittent,  or  of  decreased  accuracy  for 
users  at  these  latitudes.  The  final  constellation 
configuration  is  to  be  selected  early  in  1982.  Qtef.  2l] 

Bach  MAVSTAR  satellite  has  the  following 
characteristics:  weight  of  982  pounds,  length  of  17.5  feet 
froa  tip  to  tip  of  the  5  square  solar  panels,  three  atonic 
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frequency  standards  (clocks)  ,  and  three  nickel-cadmium 
batteries  for  operations  in  darkness  and  during  peak  load 
periods.  Stabilization  of  each  HA7STAR  is  maintained  by 
four  skewed  reaction  wheels  which  produce  three-axis 
stabilization.  MA7STAR  station-keeping  and  momentum  daaping 
capability  is  provided  by  a  hydrazine  propulsion  systea 
within  each  satellite.  Life  expectancy  of  each  unit  is  five 
years,  and  expendables  are  expected  to  last  for  up  to  7 
years  with  noraal  usage. 

Figure  3.3.  illustrates  the  parts  of  each  NA7STAR.  A 
twelve  element  shaped-beaa  helix  antenna  transaits  the 
navigation  aessage  on  two  L  band  frequencies  to  user  systea 
segments.  Simultaneously  it  receives  satellite  status, 
clock  corrections,  ionospheric  data,  and  epheaeris  constants 
frca  the  control  segment  on  two  S  band  frequencies. 
(Ref.  22]  The  two  coherent  L  band  frequencies  are  centered 
on  1227.6  MHz  (LI),  and  1575.4  MHz  (L2)  .  [Ref.  23] 
(Ref.  24]  Both  frequencies  are  necessary  to  determine  the 
total  ionospheric  delay  correction  for  signal  propagation 
time  (Table  I  ).  Zrrors  produced  by  the  ionospheric  delay 
in  signal  velocity  contribute  most  of  the  error  in  the 
pseudo-range.  By  utilizing  the  LI  and  L2  frequencies  and 
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TABLE  I.  Ionospheric  Group  Delay  as  a  Function  of  Frequency 


R  .  A  ,  B  ,  C 
C  2  1  ^ 


R  A 

C  f~2 


(neglecting  third  order  and  higher  terms  which  contribute  minimal 
error) 


where  A  =  ionospheric  delay  constant  of  conditions 
R  =  true  range 
C  =  speed  of  light 
f  =  carrier  frequency 

B  <*  (averaged  earth  magnetic  field  strength  along 
the  path)J 

=  group  delay  between  points  G  and  D 


Term  Heolected: 


Range  Error  at 
f  =  1  .5  GHz 


^  1  inch 


-  3  inches 


1 


the  following  algorithm,  the  error  is  decreased  but  not 
totally  eliminated.  [Bef.  25] 

T  '  -  T  a  A  T 

GDLi  GDL2  ' 


Where:  A 

TGDi_r  t 

L1 

where  A  =  ionospheric  constant  and  t  qdl^»  t  GDL^  are  the 
respective  ionospheric  group  delays  at  frequencies  Li  and 
L2.  Both  frequencies  are  also  nodulated  to  transmit  the  GPS 
time  codes  and  the  navigation  message. 

GPS  time  is  transmitted  using  two  pseudo-random  noise 
chip  codes — the  Coarse  Acquisition  (C/A)  and  the  Precision 
(P)  codes.  Modulation  on  the  Li  frequency  carries  both  the 
C/A  and  P  codes  in  phase  quadrature,  while  the  L2  frequency 
carries  only  the  ?  code.  These  codes  identify  each  acquired 
satellite  by  matching  the  unique  pseudo-random  noise  code 
pattern  generated  by  each  satellite  with  similar  user 


generated  codes 


Both  codes  are  required  to  measure  the  navigation  signal 
propagation  time.  The  phase  shift  is  measured  and  used  to 
aatch  the  user  generated  signal  with  the  inconing  satellite 
signal.  (Bef.  26]  The  C/A  code  has  a  repetition  period  of 
one  Billisecond  and  has  a  short  code  stream,  transmitting 
1.023  aillion  bits  per  second  (bps).  This  code  is  easily 
acquired  and  matched  with  the  user  generated  pseudo-random 
noise  code  but  only  provides  a  coarse  time  signal  for 
computing  the  pseudo-range.  However,  the  ?  code  has  a  long 
code  stream,  repeating  only  every  seven  days,  and  it  is 
transmitted  at  the  higher  rate  of  10.23  million  bps.  with 
this  higher  data  transmission  rate  and  long  repetition 
period,  a  more  accurate  time  signal  is  available  for 
pseudo-range  determination.  However,  the  long  period  and 
high  data  rate  makes  the  P  code  difficult  to  match  with  the 
user-generated  pseudo-random  noise  code  pattern  unless  the 
receiver  has  a  highly  accurate  tine  standard,  and  the 
approximate  receiver  location  is  known.  in  lieu  of  this, 
the  usual  method  of  acquiring  the  P  code  is  by  matching  the 
C/A  code  of  the  desired  NA7STAR,  then  using  the  Handover 
Word  (HOW) ,  to  locate  the  correct  P  code  sequence  within  the 
7-day  pseudo-random  noise  code.  GPS  time  is  transmitted 


with  the  HOW  every  six  seconds  and  provides  very  accurate 
GPS  tiae  to  the  user.  Both  of  these  codes  provide  the 
proper  tiae  reference  for  computation  of  the  pseudo-range. 
The  C/A  code  provides  a  less  accurate  tiae  reference  than 
the  P  code.  The  selection  of  the  code  depends  only  on  the 
navigational  accuracy  desired. 

The  navigation  message  is  also  modulated  on  frequencies 
LI  and  L2  but  at  a  auch  slower  rate  (only  50  bps)  than 
either  the  C/A  or  P  code.  On  both  frequencies,  the 
navigation  message  is  identical  and  consists  of  5  subsets  of 
6  seconds  duration  each.  2ach  navigation  message  begins 
with  the  telemetry  message  (TLtl)  and  the  HOW  message,  which 
are  used  to  transfer  from  the  C/A  to  ?  codes.  The  remainder 
of  the  entire  navigation  message  allocated  to  each 
subsection  is  transmitted  after  the  TLH  and  HOW  messages. 
Information  in  the  navigation  message  includes  approximate 
satellite  ephemerides,  status  of  other  NA7STAR  satellites, 
parameters  for  clock  corrections,  atmospheric  signal  delay 
corrections,  satellite  performance  status,  momentum  dump 
status  since  last  upload,  tiae  since  last  upload,  almanac 
and  identification  codes  for  all  satellites,  and  provisions 
for  inclusion  of  any  special  or  important  messages  from  the 
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control  system  segment  to  the  user.  [Ref.  27,  28J 
Navigation  is  thus  achievable  by  combining  the  information 
transmitted  by  the  space  system  segment  in  the  form  of  the 
C/A  and  P  codes  and  the  navigation  message. 

Accuracy  and  proper  function  of  GPS  are  governed  by  the 
second  segment — the  control  segment.  The  aontrol  segment 
monitors  and  uploads  the  individual  NAVSTAB  satellites, 
checks  operational  characteristics,  determines  the  epheaeris 
and  almanac  data  and  corrections  required,  computes  the  time 
delay  and  clock  corrections  for  each  satellite  atomic 
standard,  determines  the  atmospheric  delay  corrections  for 
pseudo-range  computations  and  includes  special  messages  in 
the  body  of  the  navigational  message.  These  functions  are 
under  the  command  and  operational  control  of  the  Master 
Control  Station,  currently  located  at  Vandenberg  AFB, 
California.  Four  unmanned  remote  monitor  stations,  located 
in  Hawaii,  Guam,  Elmendorf  AFB,  Alaska,  and  Vandenberg  AFB, 
collect  the  pseudo-ranges  to  each  satellite,  the  change  in 
pseudo-range  of  each  satellite  signal,  local  meteorological 
data,  and  the  remote  monitor  station  atonic  standard 
parameters.  Bach  set  of  <.a  is  transmitted,  on  request,  to 
the  Master  Control  Station,  where  all  the  data  is  combined 


and  the  satellite  upload  corrections  and  paraaeters  are 
computed.  These  corrections  and  paraaeters  are  then  relayed 
to  the  Opload  Station,  Vandenberg  APB,  for  transmittal  to 


the  respective  NAVSTAR  satellites.  The  uploads  occur  when 
the  satellites  are  first  visible  to  the  Opload  Station  and 
when  the  Oser  Range  Error  exceeds  4  meters,  as  computed  by 
the  flaster  Control  Station.  [Ref.  29,  30,  3lJ 

The  last  segaent,  the  user  systea  segment,  combines  the 
information  provided  to  the  space  systea  segment  (by  the 
control  segaent)  with  the  information  received  from  the 
sp^ce  system  segment  to  obtain  the  pseudo-ranges  from  four 
NAVSTAR  satellites.  These  ranges  are  computed  using  the 
following  algorithms: 

ai  =  “  cAtAi 

where  &j  =  true  slant  range  to  satellite  i,  c  *  speed  of 
light  in  a  vacuum,  tp  =  time  of  received  GPS  signal  by  user 
receiver  (assumed  to  be  simultaneous  from  all  satellites), 
t.  a  tiae  of  transmission  at  satellite  i,  and 
^tA  =  atmospheric  propagation  delay  time  from  satellite  i 
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be  expressed  in  teres  of  a  pseudo-range: 


V  58  Hi  +  cAt  a.  *  c(AVAts. * 

where  Bp  =  pseudo-range  including  clock  error,  At0  =  time 

i  K 

delay  in  the  user  receiver  cf  the  received  signal,  and 
A 1 5^.  *  tiae  delay  in  satellite  transmission  froa  satellite 
i. 

To  deteraine  the  slant  ranges: 


The  following  quantities  are  known:  x  ,  y  ,  z  obtained 
froa  navigation  epheaeris  data;  A  t.  computed  using  the 
ionospheric  delay  correction  obtained  froa  observing 
frequencies  LI  and  L2  (tropospheric  delay  neglected);  &  t  $. 
obtained  froa  the  navigation  message  as  the  satellite  clock 
error.  Thus,  there  are  four  unknowns:  x,  y,  z  (user 
position  coordinates)  and  Atp  (  user  clock  error  or  bias). 
To  solve  this  set  of  equations,  four  satellites  must  be 
observed  siaultaneously.  [Bef.  32] 
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Should  only  three  satellites  be  observable,  the  altitude 
aay  be  assumed  and  the  solutions  can  be  computed  for  the 
other  two  location  variables  and  the  user  clock  error. 
Solutions  of  this  type  (altitude  hold)  significantly  degrade 
position  accuracy. 

The  above  description  is  over-siaplif ied  since  the 
relativistic  effects  of  W&75TAB  satellite  velocity,  as 
observed  by  the  user,  are  neglected  in  terns  of  the 
satellite  tine  delay  correction  and  the  atmospheric 
propagation  delay  correction.  The  actual  user  velocity,  in 
terms  of  velocity  in  the  x ,  y,  and  z  directions,  is 
calculated  using  the  Doppler  shifts  of  the  two  frequencies 
LI  and  L2  and  the  estimated  satellite  position  taken  from 
the  ephemeris  data.  With  the  relativistic  effects  as  an 
integral  part  of  the  data  set,  solution  of  the  navigation 
position  problem  becomes  more  complex.  Phase  I  receiver 
manufacturers  have  included  them  in  the  solution  of  the 
problem  by  several  different  methods.  Methods  of  solution 
depend  on  the  usage  of  the  receiver  (whether  static,  high  or 
low  dynamic)  and  the  expected  cost/benefit  ratio  for  the 
receiver  (in  terms  of  the  total  cost,  weight  and  size). 
Current  state-of-the-art  design  technology  and  navigational 
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algorithms,  developed  since  Phase  I,  will  be  included  in  the 
Phase  II  receivers  and  will  significantly  improve  accuracy 
of  future  receivers. 

B.  BANPACK  BECEIVEH 

The  Phase  I  receiver  used  for  this  research  was  the 
Texas  Instruments  5PS  Banpack,  also  known  as  the  Banpack 
Vehicular  User  Equipment  (BVUE) .  BVUE  is  a  single  channel 
low  dynamic  receiver,  usable  for  receiver  velocities  from  0 
to  25  meters  per  second.  [Hef.  33 3  System  configuration  is 
illustrated  by  figure  3.5.  Composing  the  system  are  the 
Banpack  control  data  unit  (CDU),  antenna  and  vehicle  mount, 
and  vehicle  power  converter.  The  Banpack  is  the  receiver 
unit  (fig.  3.6.)  with  provisions  to  connect  an  external  data 
output,  through  the  BVUE  Instrumentation  System  (BIS) 
interface,  figure  3.7.  Facilities  for  an  internal 
preamplified  antenna  mount,  or  an  external  preamplified 
antenna,  and  an  attachment  for  either  an  external  power 
source  or  the  integral  battery  pack  are  also  shown  (figure 
3.7).  Together,  the  BVUE  and  the  SBP/9900  microprocessor 
data  processing  unit  have  storage  capacity  for  47,360 
sixteen  bit  words.  This  storage  capacity  allows 
interrogation  of  NAVSTAB  satellites  sequentially,  with  the 
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Fig.  3.7.  Cutaway  of  MVUE  Receiver 
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receiver  computing  a  position  fix  every  24  seconds.  The 
relativistic  velocity  effects  on  the  received  signals  are 
determined  by  aatching  the  incoming  frequencies  with  the 
expected  frequencies.  The  user  velocity  relative  to  the 
satellite  is  deterained  front  the  observed  frequency  shifts. 
The  noainal  frequency  shift  for  each  satellite  transaission 
is  coaputed  using  the  siaultaneously  transmitted  epheaeris 
data.  Total  weight  of  this  receiver  unit  is  33.5  pounds, 
and  it  has  a  volume  of  17,958  cubic  cantiaeters. 

The  second  element  of  the  MVOE,  the  CDO  (fig.  3.8.),  is 
the  input/output  device  for  the  non-autoaated  Manpack  user. 
Inputs  are  entered  by  keying  the  pressure  pads  set  in  six 
rows  of  four  pads  per  row.  To  obtain  a  navigational 
position,  the  approximate  user  latitude,  longitude,  and 
Greenwich  dean  Time  (GMT)  are  keyed.  After  the  waraup 
period  of  3  to  5  ninutes,  the  satellite  signals  are 
received,  and  the  CDO  will  display  the  actual  position. 
Additional  input  aay  include  NAYSTAH  satellite  selection, 
user  altitude,  and  other  position  datuas.  Output  options 
include  position,  GMT  tine,  and  system  error  aessages. 
Other  inputs  and  outputs  are  listed  in  the  HYOE  manual. 
£  Bef.  34] 
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Automated  users  may  elect  output  either  to  a  radio 
interface  or  through  the  MIS  interface.  For  this  research, 
a  special  interface  unit,  designed  and  constructed  at  the 
Office  of  Testing  and  Evaluation  of  Systems,  NOAA,  was  used 
to  link  the  MIS  interface  board  within  the  receiver  with  a 
Hewlett  Packard  9825T  coaputer.  With  this  interface,  user 
latitude,  longitude  and  GPS  tiae  were  logged  onto  magnetic 
tape  for  future  processing.  Pigure  3.9.  and  Table  II  are 
the  engineering  documents  used  in  construction  of  the  MVUE 
interface  unit. 

Other  MVUE  elements  are  the  antenna  and  the  antenna 
external  preamplifier  vehicle  mount.  Located  within  the 
antenna  are  two  antenna  subassemblies  associated  with 
frequencies  LI  and  L2.  Both  of  the  subassemblies  are 
connected  to  the  external  preamplifier  with  two  lengths  of 
coaxial  cable,  enclosed  in  a  flexible  conduit.  The  entire 
unit,  antenna  and  preamplifier,  is  connected  to  the  Manpack 
receiver  with  a  power  cable  and  a  coaxial  antenna  cable. 
Cable  length  is  determined  by  the  location  of  the  external 
antenna  but  is  limited  by  the  minimum  strength  required  by 
the  Manpack  receiver.  The  combined  antenna  and  preamplifier 
unit  is  36  inches  long,  and  it  weighs  1.9  pounds. 


Circuit  Diagram  of  Receiver  Interface 
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The  power  converter  must  supply  24  volts  DC  to  the 
Manpack  receiver.  During  this  research,  a  power  supply 
which  converted  120  volts  1C  to  24  volts  DC  was  used.  This 
power  supply  was  considered  more  stable  than  the  power 
converter  provided  with  the  HVtJE,  and  the  voltage  output  and 
power  were  controlled  and  aonitored  during  the  testing 
period. 

C.  DIPFEHENTIAL  GPS 

As  aentioned  earlier  in  this  chapter,  ataospheric 
propagation  delays  are  the  aajor  source  of  error  in  the 
navigational  equations  because  ionospheric  propagation  is 
poorly  understood.  The  tropospheric  propagation  delays  are 
neglected  in  the  navigation  equations,  and  the  resultant 
error  is  dependent  on  local  te a,>eratura,  pressure,  humidity, 
and  N AVSTAR  elevation  and  aziauth.  Tropospheric  corrections 
are  usually  neglected  if  the  observed  satellite  is  more  than 
5°  above  the  local  horizon.  [ Hef .  35] 

Additional  error  sources  are  receiver  and  satellite 
clock  perturbations,  lengths  of  signal  paths  from  different 
satellites,  aultipathing  delays,  satellite  epheaeris  errors, 
and  receiver  noise.  Combined,  these  factors  lead  to  errors 
of  3.6  to  6.3  aeters,  as  table  III  illustrates.  [Hef.  36] 


TABLE  1 1 1. RANGE  ERROR  BUDGET 


Uncorrected  Error  Source 


User  Eaui valent  Ranae  Error, 

- TC'i  — - - 

feet  meters 


SV  Clock  Errors  * 

Ephemeris  Errors  5.0 

Atmospheric  Delays  8.0  -  17.0 

Group  Delay  (SV  Equipment)  3.3 

Multipath  4.0  -  9.0 


1.5 

2.4  -  5.2 
1  .0 

1.2  -  2.7 


Receiver  Noise  and  Resolution 
Vehicle  Dynamics 


Root  Sum  Square  11.8  -  20.7  3.6  -  6.3 

♦Two  hours  after  update 
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One  aethod  used  to  reduce  the  aagnitude  of  this  error  is 
that  cf  differential,  or  translocated,  positioning.  To  use 
this  aethod,  one  receiver  is  positioned  over  a  known 
geographic  position.  Pseudo-ranges  to  each  satellite  are 
deterained,  and  a  differential  correction  is  obtained  for 
each  pseudo-range.  These  corrections  are  then  applied  to 
the  satellite  pseudo-range  observed  at  a  second  aobile 
receiver.  »ith  the  differential  corrections  applied,  the 
net  error  of  the  aobile  Gps  receiver  is  significantly 
reduced.  (Ref.  37] 
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17.  TEST  PROCEDPRES 


4.  POSITIONING  NETWORK 

One  of  the  considerations  in  planning  this  test  was  to 
examine  the  various  methods  of  positioning  the  shore-based 
equipment.  It  was  desirable  to  use  closely  spaced  (1-2  km) 
locations  within  the  immediate  surroundings  of  Monterey 
harbor  for  reasons  of  ready  access  by  personnel  during 
installation  of  equipment  and  for  optimum  accuracy  of  the 
positioning  systems  chosen  for  control.  4s  mentioned 
previously,  the  positioning  systems  chosen  for  the 
comparison  (the  three  theodolites  and  the  ARTEHIS)  give  best 
results  at  relatively  short  ranges. 

Two  predetermined  station  marks  were  selected  for  the 
test— 1)  USE  HON  1978,  a  bronze  disk  set  by  the  4rmy  corps 
of  Engineers,  located  on  a  sand  dune  on  the  Naval 
Postgraduate  School  property  400  meters  from  the  ocean  at  an 
elevation  of  40  meters  above  sea  level  and  2)  HOSSEL  1932,  a 
bronze  disk  set  by  the  Coast  and  Geodetic  Survey  in  a  rock 
outcrop  7  meters  above  sea  level  at  the  shoreline  on  the 
property  of  the  Hopkins  Marine  Laboratory  of  Stanford 
University.  The  third  station  mark  was  3et  in  the  center  of 


the  Monterey  coast  Guard  Pier  for  purposes  of  the  test. 
This  was  a  DMA  bronze  disk  staaped  GEOCEIVER  STA  31370. 
This  location  was  chosen  for  access  to  power  facilities 
available  at  the  pier  and  because  of  central  location  to  the 
test  area. 

The  position  of  the  Geoceiver  station  was  established 
using  a  Magnavox  Geoceiver  (Model  AN/PRR  14)  obtained  from 
the  Department  of  Satellite  Geophysics,  Satellite  Tracking 
Branch,  DMA/Hydrographic  Topographic  Center.  The  eguipaent 
was  operated  in  accordance  with  DMA  requirements  to  obtain  a 
documented  first  order  station.  The  data  were  submitted  to, 
and  processed  by,  DHA/HTC  in  Brookmont,  Maryland,  using 
precise  epheaeris  information  to  provide  the  WGS-72 
position.  Michael  Sllett  of  DMA  provided  the  technical 
expertise  for  acquisition  of  the  data  and  maintenance  of  the 
geoceiver  during  the  test  period. 

A  third  order.  Class  I,  closed  traverse  was  conducted 
following  the  standards  of  the  National  Geodetic  Survey 
(NGS)  ,  using  a  wild  T-2  theodolite  to  measure  angles  and  a 
Motorola  MRA  V  Telluroaeter  to  measure  distances.  The 
TRANSIT  position  computed  by  DMA/HTC  was  used  for  the 
initial  position  of  the  survey,  and  a  Polaris  observation 


71 


using  the  theodolite  provided  the  initial  aziauth  (fig. 
4.1.).  A  total  closure  error  in  distance  of  0.061  aeters 
and  angle  of  5.4"  was  achieved.  The  positions  deterained 
froa  the  traverse  and  used  for  the  control  network  are  as 
follows: 

OSE  MON  36°  36' 08  .424"N  1 2 1°  52*  40 . 30 1"W 
H0SSEL  36°  37*21  .882"N  1 2 1°  54 •  16 . 048"W 
GBOCBI7BB  STA  31370  36°  36*  36  .246 "N  121°  53*  29.693"W 

B.  POSITIONING  EQOIPHENT 

As  previously  described,  three  positioning  systeas  were 
used  during  the  test — the  three  theodolites,  the  GPS,  and 
the  ABTENIS.  In  addition,  the  geoceiver  was  set  to  acquire 
data  during  the  shipboard  portion  of  the  test.  During  the 
tiae  of  day  that  the  GPS  satellites  were  available  (Table 
17)  only  one  or  two  observations  per  day  of  the  saae 
TBANSIT  satellite  were  possible.  The  geoceiver  data  were 
recorded  but  were  not  used  for  this  study. 


Test  Area  with  Traverse  Sketch 


This  section  describes  the  placeaent  and  operation 
o£  the  equipment  during  the  test.  One  theodolite  was  set  on 
a  tripod  at  each  of  the  three  stations,  and  each  was 
occupied  by  an  observer  and  a  recorder.  Coaaunication 
between  stations  was  aaintained  via  portable  radios. 

Because  tests  could  be  conducted  only  between  10  PH 
and  2  AH  local  aean  tiae  (LHT)  due  to  GPS  availability 
tiaes,  it  was  necessary  to  plan  a  systeas  of  lights  for  the 
stations.  The  three  theodolite  sites  were  aade  visible  by 
placing  battery  powered  lights  directly  over  the  station 
aarks.  In  addition,  battery  operated  light  packs  for  the 
theodolites  were  used.  All  shorebased  lights  and  light 
packs  for  the  test  were  provided  by  the  Pacific  Harine 
Center,  NOAA. 

The  ship  systea  was  designed  with  an  oanidirectional 
pair  of  lights  anchored  around  the  support  pipe  for  the  GPS 
antenna.  Also,  the  equipaent  shelter  on  the  boat  deck  of 
the  ship  was  lighted  by  spotlights  so  that  the  GPS  antenna 
was  easily  seen  by  the  shorebased  theodolite  operators. 

During  the  test,  each  theodolite  observer  initialled 
on  one  of  the  other  stations  and  recorded  the  reading  on  a 
log  sheet.  The  theodolite  operators  then  intersected  the 
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lights  on  the  ship  GPS  antenna  using  the  theodolite 
telescope  and  tracked  the  ship  until  a  'nark*  was  called 
over  the  radio  fros  the  personnel  aboard  the  ship.  The 
angle  measured  at  this  tiae  was  recorded  on  the  log  sheet  as 
a  fix.  The  tiae  interval  between  fixes  was  approximately 
one  ainute.  Each  series  of  approximately  twenty  fixes 
constituted  a  set.  At  the  end  of  each  set,  a  final  reading 
was  aade  of  the  angle  to  the  reference  station.  The  fixes 
were  nuabered  consecutively  throughout  the  test,  and  the 
Greenwich  Hean  Tiae  (GMT)  of  each  fix  was  recorded  on  the 
ship  log  sheets.  The  length  of  the  sets  depended  on  the 
ship  position  and  heading. 

To  correct  the  observed  theodolite  angles,  the  mean 
of  the  final  pointing  and  the  initial  pointing  was 
determined  for  each  set  and  for  each  station.  This  value 
was  assuaed  to  be  a  constant  throughout  the  set  and  was 
reaoved  from  each  angle.  A  colliaation  error  for  each 
theodolite  was  also  reaoved  froa  each  theodolite 
observation.  The  sane  observers  and  theodolites  were  used 
nightly  at  all  stations  except  MOSSEL,  where  two  different 
individuals  observed  over  the  test  week.  This  continuity  of 
personnel  and  eguipaent  was  designed  to  reduce  operator  and 
instruaent  error  to  a  minimum. 
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2.  2££  aaa  ARTS  MIS 

To  colocate  the  different  pieces  of  equipment  at  the 
GEOCSIVER  STA  31370,  a  wooden  tower  was  built  as  shown  in 
figure  4.2.  This  tower  was  set  over  the  station  for  the 
period  of  the  test.  An  UPS  equipment  vehicle  was  parked  as 
close  as  possible  to  the  tower.  Electrical  power  was 
obtained  from  the  Honterey  coast  Guard  facility. 

Prior  to  the  ship  test,  both  GPS  antennas  were 
placed  on  the  tower  and  centered  over  the  station  mark.  Por 
the  first  two  nights,  9  and  10  Ray,  equipment  problems 
allowed  the  use  of  only  one  system  each  night  and  data  were 
acquired  only  to  determine  the  stability  of  each  system. 
During  the  next  night,  11  Hay,  position  outputs  for  the  two 
simultaneously  operating  GPS  receivers  were  recorded.  This 
static  test  gave  a  determination  of  any  difference  between 
the  two  receivers  under  almost  identical  conditions. 

Por  the  ship  test,  the  GPS  antenna  was  centered  over 
the  mark  and  was  anchored  at  the  top  of  the  tower;  the 
geoceiver  antenna  was  centered  over  the  mark  and  anchored  at 
the  middle  level;  and  the  theodolite  on  the  tripod  was  set 
over  the  mark  at  ground  level.  This  arrangement  provided 
colocation  of  all  three  systems.  It  was  assumed  that  no 
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Antenna  Tower  on  Coast  Guard  Pier 


error  was  introduced  by  the  presence  of  the  wooden 
structure.  The  sho rebased  ABTEHIS  antenna  unit  was  offset 
frca  the  theodolite  to  eliminate  interference  in  the  signal 
due  to  the  tower.  An  offset  position  was  computed  relative 
to  the  station  mark  and  was  used  in  postprocessing  of  the 
ARTEMIS  data. 

Except  for  the  antennas  for  the  GPS,  the  ARTEMIS, 
and  the  geoceiver  systems,  all  other  eguipaent  was  installed 
in  the  eguipaent  vehicle,  where  an  operator  was  stationed 
throughout  the  test  to  monitor  the  data  quality. 

The  GPS  Manpack  is  a  manually  operated  system.  As 
mentioned  previously,  it  was  decided  that  manual  operation 
of  GPS  was  detrimental  to  the  goal  of  the  test,  and  that  an 
interface  that  would  provide  an  automated  printout  of 
positions  from  the  Manpack  as  well  as  a  printout  of  GPS  time 
would  be  desirable.  Knute  Berstis  and  Gary  Whitsell  of  the 
Office  of  Testing  and  Evaluation  of  Systems  (OTES) ,  NOAA, 
assembled  two  interfaces  that  provided  data  storage  on 
cassette  tapes  as  well  as  the  desired  printout.  Each  systea 
operated  in  conjunction  with  a  HP-9825  computer  and  a 
HP-7845  printer/plotter.  Software  was  provided  by  OTES, 
NOAA.  This  arrangement  permitted  storage  of  the  satellite 


messages  for  further  processing  and  provided  an  autonated 
printout  of  the  ARTEMIS  data.  A  GPS  Manpack  had  been  sent 
to  OTES  for  testing  during  the  construction  of  the 
interface.  Technical  expertise  and  operation  of  the 
interface  during  the  test  were  provided  by  OTES. 

The  data  from  the  GPS  receiver  were  printed  at  an 
interval  of  24  seconds,  and  consecutive  data  were 
continually  recorded  on  tape  and  printed  during  the  test 
periods.  In  addition,  a  annual  record  of  the  GPS  position 
froa  the  COO  w as  logged  as  the  • marks*  were  called  via  the 
radios. 

The  ship  positioning  systems  were  installed  as  shown 
in  figures  4.3.  and  4.4.  An  equipment  shelter  was  secured 
to  the  boat  deck  of  the  R/V  ACAMIA.  A  length  of  six-inch 
diameter  polyvinyl  chloride  (PVC)  pipe  was  attached  to  the 
aft  wall  of  the  building  on  the  centerline  of  the  ship.  The 
GPS  antenna  and  the  system  of  lights  previously  mentioned 
were  set  on  top  of  the  pipe. 

The  ARTEMIS  mobile  antenna  unit  was  set  on  a 
platform  located  on  the  centerline  of  the  ship  and  at  the 
top  of  the  building  high  enough  to  not  interfere  with  the 
signal  froa  the  ship  radar  beam.  Por  the  period  of  the 
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test,  the  weather  conditions  allowed  ship  operations  without 
use  of  the  radar,  thus  there  was  no  possibility  of  damage  to 
the  AHTEMIS  unit. 

All  other  equipment — the  AHTEMIS  receiver  unit,  the 
GPS  receiver,  control  data  unit,  power  filter,  and 
interface,  the  HP-9825,  the  HP-7845  plotter/  printer,  and 
the  real  time  clock — was  placed  inside  the  building.  An 
operator  monitored  the  equipment  during  the  ship  tests  and 
another  individual  radioed  the  'mark*,  or  fix,  to  the 
theodolite  operators  on  shore.  The  equipment  used  at  the 
shore  location  was  installed  in  a  vehicle  and  consisted  of 
the  same  components  as  on  the  ship.  Operation  of  equipment 
duri..  j  the  ship  tests  was  monitored  by  an  individual  located 
in  the  vehicle. 

C.  LOG  OP  EVENTS 

The  test  was  divided  into  two  portions — 1)  a  shorebased 
test  which  consisted  of  two  nights  of  stability  testing,  and 
one  night  of  receiver  comparison  tests,  and  2)  the  ship  test 
consisting  of  three  nights  of  shipboard  operation.  The  H/V 
ACANIA  from  the  Naval  Postgraduate  School  was  the  platform 
used  for  the  shipboard  tests. 
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1.  Stability  Jafii 

Operations  began  on  9  May,  1981,  GMT.  One  of  the 
GPS  systems  was  placed  over  the  station  mark,  and  continuous 
data  was  recorded  during  the  satellite  availability  periods. 
In  this  manner,  system  1  was  tested  on  9  May,  and  system  2 
was  tested  on  10  May.  A  position  was  defined  as  a  quality 
observation  if  four  or  more  satellites  were  used  by  the 
receiver  to  obtain  the  output. 

2.  Static  Test 

Both  GPS  antennas  were  centered  over  GEOCEIVEB  ST fc 
31370,  and  data  from  simultaneous  operation  of  both 
receivers  were  recorded  at  the  equipment  vehicle.  One 
antenna  was  placed  on  top  of  the  wooden  structure  and  the 
other  immediately  below  it  on  the  middle  level.  The 
vertical  displacement  was  assumed  to  cause  negligible  error 
in  the  horizontal  positions.  Data  ware  acquired  on  11  May 
free  both  receivers  with  a  quality  observation  defined  as  a 
position  determined  from  four  or  more  satellites. 

3.  shipboard  Test 

For  the  entire  period  of  shipboard  testing,  the  wind 
speeds  were  less  than  5  knots  and  the  seas  were  calm  with  1- 
2  feet  or  less  of  swell.  Visibility  was  excellent  throughout 
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the  test  week.  The  ship  track  for  each  set  is  shown  in 
figure  4.5.  A  log  of  subsequent  operations  follows: 

12  May - The  GPS  Manpack  receiver  was  installed  on  the  R/V 

ACANIA  and  GPS  and  theodolite  operations  were  conducted. 
The  ship  was  allowed  to  drift  during  sets  1-3  (fix  1-69) 
with  23  fixes  per  set,  in  a  east-northeasterly  direction. 

13  May - The  ship  cruised  at  4  knots  on  track  lines  froe 

south  to  north  for  sets  4  and  6,  and  lines  from  north  to 
south  for  sets  5  and  7.  During  fixes  139  through  157  the 
radio  at  the  Coast  Guard  pier  was  disabled,  and  no  data  were 
acquired  by  either  the  theodolite  observer  or  the  operator 
in  the  equipnent  vehicle. 

14  May - The  shipboard  GPS  system  was  inoperable.  The 

trouble  was  traced  to  a  faulty  antenna  cable,  requiring 
replacement.  Mo  data  were  acquired. 

15  May - The  ARTEMI S,  the  GPS  and  the  theodolite  network 

were  all  used  for  this  test  period.  Data  from  sets  9 
through  13  were  acquired  (fixes  165-262)  . 

16  May - The  control  data  unit  on  the  shore  based  GPS  system 

failed  before  the  start  of  ship  operation  and  could  not  be 
repaired  without  factory  maintenance. 
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A  total  of  252  fixes  were  recorded  by  the  theodolite 
network  and  all  GPS  positions  corresponding  to  the 
theodolite  fixes  were  aanaally  logged.  The  ARTSHIS  aziauth 
and  distance  inforaation  and  consecutive  GPS  positions  on 
the  ship  and  in  the  vehicle  were  stored  on  cassette  tapes 
and  printed  by  the  HP-7845  printer/plotter. 
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7.  DATA  PROCESSING  METHODS 

Processing  of  the  thesis  data  was  divided  into  four 
interrelated,  but  separate,  tasks.  These  tasks  were  the 
basic  geodetic  survey,  theodolite  and  ARTEMIS  ship  position 
deteraination,  differential  SPS  ship  position  deter ainat ion, 
and  the  statistical  analysis  of  the  ship  position 
ccaparisons.  Each  task  will  be  discussed  separately. 

A.  BASIC  SUR7EY  COMPUTATIONS 

As  described  earlier  in  the  Test  Procedures,  a  basic 
geodetic  survey  of  the  test  area  was  conducted  to  place  the 
theodolite  control  stations  on  the  Horld  Geodetic  System  of 
1972,  or  WGS-72.  To  begin  the  survey,  a  Magnavox  AN/PER  14 
Geoceiver  was  placed  over  the  center  survey  disk,  GEOCEIVER 
STA  31370  1981,  and  the  digital  Doppler  data,  from  TRANSIT 
satellite  68,  were  recorded  onto  paper  tape.  Thirty-three 
passes,  over  a  period  of  two  and  one  half  weeks,  were 
recorded,  but  analysis  of  the  equipment  maintenance 
diagnostics  indicated  a  possible  error  in  the  records. 
These  data  were  therefore  discarded. 

A  second  set  of  forty  satellite  passes  of  TRANSIT 
satellite  68  were  recorded  over  a  ten  day  period  and  were 
determined  to  be  of  acceptable  quality  according  to  the 
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standards  of  the  DMA  Department  of  Geodesy.  This  second  set 
of  data  was  then  submitted  to  DMA ,  Department  of  Geodesy, 
Computation  Section,  for  analysis.  Also  submitted  were 
meteorological  data  taken  during  the  passes  and  additional 
information  required  for  documentation  of  a  first  order 
IHAHSIT  position.  Satellite  68  Doppler  data  for  the  survey 
was  then  processed  in  conjunction  with  the  tracked  satellite 
position  by  DMA  computer  facilities.  A  final  geodetic 
position  on  the  HGS-72  ellipsoid  was  obtained  for  GEOCEIVER 
STA  31370.  Coordinates  of  the  station  are  listed  in  the 
Test  Procedures  section.  The  computed  height  relative  to 
the  ellipsoidal  surface  is  *33.39  meters.  Computation  of 
the  geographic  position  utilized  a  minimum  of  38  individual 
THAMSIT  satellite  passes  recorded  for  GEOCEIVEH  STA  31370. 
This  provided  for  an  accuracy  of  one  part  in  10s  (first 
order  accuracy  standard).  [Bef.  38] 

To  orient  the  positioning  network,  an  initial  azimuth 
from  GEOCEIVER  STA  31370  to  station  MUSSEL  1932  was  computed 
from  two  sets  of  observations  of  POLABIS.  All  observations 
and  computations  of  the  star  sighting  utilized  procedures 
outlined  in  the  national  Ocean  Survey  Photogrammetric 
Instructions  Humber  4,  Revision  1,  4/24/73,  the  Hydrographic 
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Manual,  and  the  Manual  of  Geodetic  Triangulation,  SP-247. 

[Ref.  39,  403 

Having  the  initial  azimuth,  distances  vere  aeasured 
between  all  three  control  stations  using  a  Telluroaeter  NBA 
5  aicrovave  distance  measuring  unit.  Twenty  observations  of 
distance  were  made  for  each  of  the  three  lines  aeasured: 
GEQCEIVER  STA  313 70  to  MUSSEL,  G30CEIVBR  STA  31370  to  OSE 
MON,  and  MOSSEL  to  OSE  MON.  The  average  of  the  distance  for 
each  line  was  then  corrected  for  instrument  error  and 
meteorological  conditions  following  the  instructions 
contained  in  the  Operators  Manual  and  the  Telluroaeter 
Manual.  [Ref.  41]  Resulting  slope  distances  vere  then 
reduced  to  geodetic  distances  using  the  respective 
corrections  for  slope  to  horizontal  distance,  reduction  of 
elevation  to  sea  level,  and  the  chord  to  arc  distance. 
These  corrections  were  computed  using  the  WG S  1972 
ellipsoidal  parameters  and  the  methods  outlined  in  the 
Hydrographic  Manual,  the  Manual  of  Geodetic  Triangulation, 
and  Introduction  to  Surveying.  [Ref.  42]  The  resulting 
geodetic  distances  meet  third- order,  class  I  specifications 
for  traverse  accuracy,  1:20,000.  In  addition,  angles  to  the 
other  stations  were  observed  at  each  control  station,  giving 
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a  total  o£  three  sets  of  angles.  Each  set  vas  coaprised  of 
four  direct  and  four  reverse  plate  observations  with  a  wild 
T-2  theodolite.  Each  pair  of  direct  and  reverse 
observations  was  corrected  for  initial  pointing  errors. 
Observations  and  computations  were  Bade  using  the  methods 
detailed  in  the  Hydrographic  Manual,  the  Manual  of  Geodetic 
Triangulation,  and  the  Manual  of  second  and  Third-Order 
Triangulation  and  Traverse.  [Bef.  43]  The  resulting 
closure  meets  the  specifications  for  traverse  accuracy  of 
three  seconds  per  traverse  station  or  a  closure  of  6 
seconds. 

With  the  geodetic  position  of  GEOCEIVER  STA  31370  known, 
as  well  as  the  azimuth  from  GEOCEIVER  STA  31370  to  MUSSEL 
and  the  distance  from  GEOCEIVER  STA  31370  to  MUSSEL,  the 
geographic  position  of  MUSSEL  was  computed  on  the  WGS  1972 
ellipsoid  using  the  direct  geodetic  computation.  This 
computation  used  the  Rainsford’s  method  modified  by  T. 
Vincenty  and  including  the  Helmert's  elliptical  terms. 
[Ref.  44] 

Next,  using  the  azimuth  from  MUSSEL  to  USE  MON,  the 
geodetic  distance,  and  the  results  of  the  direct  geodetic 
computations,  the  geographic  position  of  USE  MON  was 


determined.  Also,  the  ARTEMIS  geographic  position  was 
calculated  as  an  offset  using  an  aziauth  and  distance  fron 
GE0C3IVEH  STA  31370.  With  these  coaputations,  the  geodetic 
survey  was  coapleted. 

Positioning  of  the  ship,  3/V  ACANIA,  required  one 
aziauth,  or  line-  of-position,  froa  each  of  the  control 
stations  MUSSEL,  GEOCEIVER  STA  31370,  and  USE  SON.  These 
lines  of  position  were  used  to  triangulate  the  vessel 
position.  Angles  froa  each  of  the  three  theodolite 
locations  were  observed  as  described  in  the  Test  Procedures. 
Each  of  the  three  angular  aeasureaents  were  then  corrected 
for  mean  pointing  errors  and  the  individual  theodolite 
colliaation  errors.  Each  corrected  angle  was  then  adjusted, 
by  control  station,  to  obtain  an  aziauth  froa  south,  as 
required  to  triangulate  the  vessel  position.  The  three 
resulting  aziauths  and  the  corresponding  control  station 
latitude  and  longitude  were  then  used  in  a  least  squares 
three  line  of  position  intersection  computation.  Results  of 
the  coaputation  were  a  vessel  position  in  latitude  and 
longitude,  error  ellipse  seai-aajor  radius,  error  ellipse 
seai-ainor  radius,  aziauth  froa  south  of  seai-aajor  radius, 
and  circular  standard  error.  [Ref.  *5]  Figure  5.1 
illustrates  these  error  ellipse  parameters. 
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THE  ERROR  ELLIPSE 


The  mean  displacement  of  a  point  to  be  expected  in  the  direction  of 
the  two  coordinate  axes,  X  and  V,  due  to  measuring  errors,  can  be  ex¬ 
pressed  by  the  quantities  and  cr 

When  the  north-south  axes  are  rotated,  two  axes  (X  and  Y)  perpendicular 
to  one  another  are  obtained,  corresponding  to  the  maximum  and  minimum 
values  of  cr  and  cr  The  maximum  and  minimum  values  of  o'  ,  and  the 
clockwise  rotation  angle  (o<)  of  the  major  axes  of  the  error  ellipse  are 
included  as  part  of  the  output  of  the  Geodetic  Package  and  can  be  used 
as  input  if  desired. 


The  circular  standard  error  (CSE)  is  also  included  in  the  output  with  the 
error  ellipse  parameters.  The  value  of  CES  is  the  radius  of  a  circle 
whose  area  is  equal  to  the  area  of  the  error  ellipse.  That  is: 


f  fP 


=  or 

max  min 


Fig.  5.1  . 
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For  this  study,  262  vessel  positions,  or  fixes,  were 
observed  and  recorded.  On  exaeination  during  processing,  51 
fixes  were  found  to  contain  errors  and  data  oiissions  which 
rendered  the  vessel  positioning  computation  invalid  or 
indeterminate.  The  remaining  211  fixes  comprised  the  data 
set  which  was  then  broken  down  into  smaller  sets  for 
analysis  and  comparison.  These  sets  were  numbered  1  through 
13  and  categorized  according  to  vessel  track.  All  sets  were 
then  processed  using  several  different  methods  to  determine 
the  most  accurate  approach. 

B.  THEODOLITE  AND  ARTEMIS  SHIP  POSITION  DETERMINATION 

Before  processing  the  OPS  Manpack  data,  the  ARTEMIS  data 
were  processed  to  ensure  the  vessel  positions,  as  obtained 
by  theodolite,  were  correlated  with  the  correct  automated 
and  hand  recorded  S PS  data.  ARTEMIS  data  consisted  of  an 
angle,  MUSSEL-ARTEMIS  OFFSET- vessel,  and  distance  from  the 
ARTEMIS  OFFSET  to  the  vessel.  The  angular  measurement  was 
first  converted  to  an  azimuth  from  south  using  the  azimuth 
of  ARTEMIS  OFFSET  to  MUSSEL  as  the  correction  and  adding  the 
measured  angle.  Vessel  position  based  on  ARTEMIS  was  then 
computed  using  the  direct  geodetic  program  (also  used  for 
the  basic  survey) .  Comparisons  of  the  resulting  geographic 


positions  and  the  theodolite  derived  geographic  positions, 
were  made  for  corresponding  fix  numbers  using  the  inverse 
geodetic  computation  on  the  WGS-7  2  ellipsoid.  The  inverse 
coaputation  utilizes  two  geographic  positions  as  inputs  for 
calculating  the  forward  and  back  aziauths  and  the  distance 
between  the  positions.  [Hef.  46]  For  comparison,  if  the 
two  positions  differed  by  more  than  20  aeters,  the  ARTEMIS 
positions  on  either  side  of  the  suspect  ARTEMIS  fix  were 
then  used  to  compute  the  inverse  coaputation  to  the 
theodolite  ship  position.  In  seven  fixes,  the  fix  position 
values  were  changed  in  the  May  15  GPS  data  to  give  more 
accurate  results. 

C.  DIFFER ENTIA1  GPS  SHIP  POSITION  DETERMINATION 

With  the  check  and  adjustment  of  the  GPS  data  and  the 
theodolite  determined  vessel  position  complete,  the  first 
GPS  data  processing  was  the  comparison  of  the  uncorrected, 
non-differential  ship  GPS  Manpack  receiver  data  to  the 
theodolite  computed  ship  position.  Using  the  inverse 
geodetic  computation,  the  forward  azimuth  and  the  distance 
between  the  theodolite  computed  ship  position  and  the 
ncn-dif ferential  GPS  ship  position  were  completed.  These 
results  were  then  used  in  the  analysis  of  the  differential 
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processing  to  exaaine  the  stability  of  the  shipboard  Manpack 
receiver  compared  to  the  GPS  Manpack  receiver  located  over 
GE0CEI7E8  ST&  31370. 

There  were  two  aethods  of  analysis  of  the  differential 
GPS.  The  first  of  these  was  the  determination  of  a  static 
correction  to  be  applied  to  all  the  GPS  data.  This  static 
correction  in  the  differential  GPS  computations  was  to 
account  for  a  suspected  bias  between  the  two  similar,  but 
not  identical,  Manpack  receivers.  A  aean  static  correction 
was  computed  for  the  most  stable  reception  period  during  the 
static  tests  on  11  May.  aith  the  two  receivers  colocated 
over  the  GEOCEIVES  STA  31370,  a  mean  correction  of  -0.954 
seconds  of  latitude  and  -0.307  seconds  of  longitude  was 
found  between  GPS  Manpack  system  1  and  System  2  during  times 
07:16:00  GMT  and  07:25:12  GMT.  This  static  correction  was 
then  added  to  the  latitude  and  longitude  of  GPS  Manpack 
System  1  (located  over  GE0CEI7ER  STA  31370)  data. 

The  second  method  of  analysis  or  processing  in  the 
differential  mode  used  the  latitude  and  longitude  obtained 
directly  froa  the  GPS  Manpack  receivers.  In  either 
processing  method,  the  computation  of  the  differential 
corrector  to  the  shipboard  GPS  Manpack  receiver  data  was  the 
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saae.  An  explanation  of  the  computations  and  processing  of 
the  differential  correction  for  both  aethods  follows. 

cooputation  of  the  differential  corrector  for  each  GPS 
fix  position  was  accoaplished  by  using  the  geographic 
position  obtained  with  the  stationary  Manpack  receiver, 
located  over  G20CEIVBR  STA  31370,  and  the  geographic 
position  of  GEOCEIVER  STA  31370.  These  two  positions  were 
then  used  in  the  inverse  geodetic  coaputation  to  obtain  a 
differential  correction  vector  (forward  aziauth,  and 
distance)  froa  the  stationary  GPS  receiver  to  GEOCEI7ER  STA 
31370.  In  figure  5.2,  the  differential  correction  vector  is 
labelled  MA".  Using  the  direct  geodetic  coaputation,  the 
differential  correction  vector  was  applied  to  the 
corresponding  shipboard  GPS  Hanpack  receiver  position,  and 
the  differential  GPS  position  was  obtained.  These  esulting 
positions  were  coaputed  on  a  one-to-one  basis  using  the  fir 
nuabers  as  identifiers.  Vector  "3",  the  differential 
corrector  vector  which  is  identical  to  the  vector  "A"  in 
figure  5.2,  is  applied  to  the  shipboard  GPS  Manpack  receiver 
and  results  in  the  differential  GPS  ship  position. 

Coaparison  of  the  differential  GPS  ship  position  with 
the  theodolite-coaputed  ship  position  was  accoaplished  using 
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the  inverse  geodetic  coaputation  for  both  methods  of 
differential  GPS  processing.  The  resulting  forward  azimuth 
and  distances  for  both  processing  methods  were  used  in  the 
statistical  analysis. 

D.  STATISTICAL  DATA  AHALTSIS 

Statistical  analysis  was  done  on  data  separated  into 
sets  according  to  the  method  of  data  processing.  The  static 
data  and  stability  data  were  graphed  as  functions  of 
latitude  and  longitude  versus  time.  User  range  error, 
furnished  by  the  Space  Division,  GPS/OL-AQ,  Vandenburg,  AFB, 
was  plotted  as  a  function  of  time  for  each  day  of  the  test 
(see  Results).  Other  graphed  relationships  are  the 
gecgraphic  positions  of  the  stationary  GPS  flanpack  receiver 
data  versus  the  GE0CEI7ER  STA  31370  and  the  theodolite 
semi- major  radius  as  a  function  of  fix  number  (see  Results). 

Other  statistical  data  analyses  were  represented  by 
histograms.  Plotted  on  these  histograms  are  the  relative 
frequency  on  the  left  margin,  variable  bar  increments  scaled 
on  the  bottom,  number  of  variable  observations  in  each  bar 
increment  on  the  top,  and  the  total  sample  size  on  the  top 
margin.  Rithin  the  histogram  the  frequency  bars  are  denoted 
by  asterisks,  the  mean  by  a  capital  M ,  the  quartiles  by 
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periods  and  the  eapirical  density  function  by  the  capital  F 
curve.  The  eapirical  density  function  was  coaputed  using: 

1  v  (x.-z) 

F(Z)  _  -  r  »  —3 - 

N(B(H)  )  i=1  B  (H) 

[Bef.  47]  where:  P(z)  =  the  eapirical  density  function,  z 

*  the  freguency  of  observation,  x  *  individual  observation, 
B  *  the  total  nuaber  of  observations,  9  *  weight  function,  B 

*  bandwidth, 

i r  »  i  -  z 

B (N)  *  range  of  all  observations  /  "fiT 
The  statistical  values  listed  below  the  histograa  are 
defined  in  Table  7.  Coaputed  histograas  show  the  distances 
between:  1)  the  theodolite- coaputed  ship  position  and  the 

differential  GPS  ship  position,  both  with  and  without  the 
static  correction,  by  set;  2)  the  theodolite-coaputed  ship 
position  and  the  differential  GPS  ship  position,  without 
static  correction,  for  five  satellite  observation  periods  on 
each  of  the  three  ship  test  days;  3)  the  theodolite-coaputed 
ship  position  and  the  non-differential  GPS  ship  position  for 
the  five  satellite-observation  periods  on  each  of  the  three 


TABLE  V.  STATISTICAL  DEFINITIONS 


The  definitions  of  the  basic  statistics  computed  by  the  program  are 
listed  below:. 

MEAN:  Average  of  the  sample 

MEDIAN:  Mid-Value  of  the  sample,  if  there  are  an  odd  number  of 

sample  points,  or  the  average  of  the  two  middle  values 
for  an  even  number  of  points 

TRIMEAN:  0.25  *  (Q]  +  2Q2  +  Q3),  where  Qi  are  the  quartiles 

MIDMEAN:  The  average  of  all  sample  values  lying  between  the  upper 

and  lower  quartiles 

MIDRANGE:  Average  of  the  maximum  and  minimum 

GEOMETRIC  MEAN:  GM  .  .  x"n 

HARMONIC  MEAN:  HM  =  n  “l/x. 

VARIANCE:  Unbiased  estimators  for  variance  and  standard  deviation 

are  used  -  the  square  of  the  standard  deviation 

STANDARD  DEVIATION:  S  =  |JTx?  -  n  1ft 

COEFFICIENT  OF  VARIATION:  Standard  deviation  /  (Meanl  when  the 

mean  is  less  than  10.  **  -30,  the  coeffic¬ 
ient  of  variation  is  set  to  zero. 

MEAN  DEVIATION:  The  average  of  the  sum  of  the  absolute  differences 
between  the  sample  values  and  the  median 

RANGE:  Maximum  -  minimum 

MIDSPREAD:  Q3  -  ,  also  called  the  interquartile  distance 

M3:  Third  central  moment.  Unbiased  estimator  is  used 

M4:  Fourth  central  moment.  Unbiased  estimator  is  used 

COEFFICIENT  OF  SKEWNESS:  M3  /  (Standard  Deviation)  **  3 

COEFFICIENT  OF  KURTOSIS:  (M4  /  (Standard  Deviation)  **  4)  -  3 

BETA1 :  Biased  estimate  of  third  central  moment;  can  be  used  in 

testing  for  normality 

BETA2:  Biased  estimate  of  fourth  central  moment 

MAXIMUM:  Largest  sample  value 

MINIMUM:  Smallest  sample  value 

QUANTILES:  The  alpha-quantile,  x(alpha),  is  the  solution  to  the 

equation: 

Probability  (x  ,le.  x(alpha)  )  =  alpha 


ship  test  days;  and  4)  the  GE0CEI7EB  STA  31370  position  and 
the  uncorrected  stationary  GPS  position  for  the  three  days 
of  ship  -ests.  In  addition ,  the  aziauths  between  the 
theodolite-coaputed  static  corrections,  for  each  subset, 
were  graphed  as  a  series  of  histcgraas.  To  illustrate  the 
difference  between  the  GE0CEI7EB  STA  31370  position  and  the 
stationary  GPS  position,  the  azieuths  between  the  two 
positions  are  graphed  on  histograas  for  each  of  the  three 
test  days. 
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VI.  E8B0B  ANALYSIS 


The  error  sources  contributed  by  each  positioning  systee 
are  described.  Estimates  of  the  effect  of  each  on  the  final 
results  are  determined  where  possible. 

A.  AHTEHIS  POSITIONING  SYSTEH 

The  sources  of  error  in  the  AHTEHIS  ship  position  can  be 
divided  into  two  categories — those  sources  of  error  inherent 
in  all  aicrowave  positioning  systems  and  those  unique  to  the 
AHTEHIS  systen.  All  are  investigated  and  the  contribution 
of  each  to  the  total  error  is  dete rained, 
t.  fllstqwayg  Svstea  gCS8Cg 

Sources  of  systea  error  include  range  holes, 
inaccuracies  in  calibration,  non-ideal  geometry  of  the  lines 
of  position  describing  the  fix,  and  interference  with 
reflections  of  the  signal  froa  nearby  objects. 

Bange  holes  are  a  cause  of  unreliability  in 
aicrowave  positioning  systeas.  The  problem  exists  when  the 
sea  surface  is  cala  enough  to  reflect  part  of  the 
transaitted  signal.  The  reflected  signal  will  arrive  at  the 
receiving  antenna  at  a  later  tiae  than  the  signal  which 
travels  directly,  causing  a  phase  difference  between  the  two 
signals.  A  zone  or  region  of  destructive  interference  of 
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the  directly  transmitted  signal  and  the  reflected  signal 
vill  occur  vithin  which  the  vessel  will  receive  an 
attenuated  signal.  Rhen  the  phase  difference  of  the  two 
paths  is  exactly  180  °,  no  signal  will  be  received  provided 
the  reflected  signal  is  as  strong  as  the  direct-path  signal. 
Prediction  of  the  range  holes,  therefore,  is  useful  in 
survey  planning  so  as  to  avoid  this  problea. 

The  range  hole  is  a  function  of  shore  unit  and  ship 
unit  antenna  heights,  the  wavelength  of  the  systee,  and  the 
distance  froa  the  shorebased  unit.  If  the  operating  area  is 
known,  antenna  heights  can  be  deternined  and  equipment  can 
be  installed  so  that  no  range  holes  exist  in  the  survey 
area.  If  fixed  antenna  heights  are  known,  the  region  of 
uncertainty  can  be  deterained  froa  the  equation:  [Ref.  48] 
The  following  equation  is  used  for  this  coaputation: 

Dn  =  2(h1)  (h2)/nX, 

where  \  is  the  wavelength  (constant  for  AETBHIS  at  .032  a) , 
hi  and  h2  are  antenna  heights,  and  Dp  is  the  distance  froa 
the  shore  based  unit  to  the  nth  range  hole.  For  this  study, 
the  antenna  height  of  the  fixed  unit  above  Bean  sea  level  at 
the  GB0C3I73R  STA  31370  was  3.8  aeters.  The  aobile  antenna 
was  installed  on  top  of  the  eguipaent  shelter  on  the  boat 


deck  of  the  H/7  ACAHIA  at  a  height  of  7.9  aeters  above  the 
the  waterline.  Two  range  holes  within  the  survey  area  were 
deterained  to  be  located  at  ranges  of  1876  aeters  for  n  *  1 
and  933  aeters  at  n  *  2  (fig.  6.1).  These  two  ranges  were 
coaputed  for  the  case  of  aean  lower  low  tide.  The  change  in 
width  of  the  range  hole  due  to  a  tidal  fluctuation  of  1.5 
aeters  during  the  test  is  shown  in  figure  6.1.  (During  the 
night  of  15  Hay  the  sea  surface  was  sufficiently  disturbed 
by  a  two-foot  swell  and  a  surface  wind  of  2-4  knots  such 
that  surface  reflection  was  ainiaal  and  no  loss  of  signal 
was  observed  at  the  coaputed  ranges) . 

Calibration  of  the  ABTEHXS  systea  was  accomplished 
in  two  stages.  As  part  of  the  study,  the  ASTEHIS  was 
obtained  for  a  one  week  test  of  the  systea  froa  13  to  17 
April.  During  this  tine  it  was  calibrated  for  aziauth  and 
range  over  the  baseline  between  the  SEOCEIVEB  STA  31370  and 
HOSSBL.  The  aziauth  readings  reguired  no  adjustment  but  the 
range  was  reduced  by  two  aeters.  The  eguipaent  was  returned 
to  Marinav  Corporation  to  fulfill  prior  coaaitaents  and  then 
was  returned  to  Honterey  for  the  ship  test  on  15  Hay. 
Before  the  evening  operation,  the  ABTEHIS  was  used  to 
aeasure  the  saae  baseline  as  previously.  A  plus  two  aeter 
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Fig. 6.1.  ARTEMIS  Range 


bias  was  observed  at  this  time,  bat  no  immediate 
coapensation  was  aade.  This  error  was  assumed  constant 
throughout  the  test  and  was  included  in  the  calculation  of 
che  ship-determined  ARTEMIS  position  during,  postprocessing. 

One  of  the  advantages  of  the  ARTEMIS  system,  unlike 
other  microwave  positioning  systems,  is  that  the  geometry  is 
always  ideal  for  the  lines  of  position  that  determine  the 
fix.  The  range-azimuth  technique  of  determining  the  vessel 
position  defines  an  uncertainty  ellipse  having  axes  that  are 
always  mutually  perpendicular.  This  allows  prediction  of 
errors  as  a  function  of  range  only,  independent  of  the 
angle. 

The  manufacturer's  specified  error  in  range 
measurement  of  the  ARTEMIS  is  +  1.5  meters  for  all  line  of 
sight  ranges.  One  axis  of  the  uncertainty  ellipse  is 
therefore  constant.  The  other  axis,  that  describes  the 
error  introduced  by  the  system  azimuth  uncertainty,  is  60  cm 
per  1000  a  range.  (This  is  determined  from  an  angular  error 
of  0.033  °  or  2*  of  arc).  Both  system  error  values  represent 
an  accuracy  of  twice  the  standard  deviation  value,  or  a  95< 
probability  of  the  true  value  occurring  within  the  margin 
allowed  by  these  parameters. 
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The  positional  uncertainty  can  be  expressed  in  teres 
of  the  radius  of  a  circle  of  equivalent  probability  within 
which  the  true  fix  is  expected  to  occur.  The  eethod  of 
coaputation  of  this  radius  is  described  in  Appendix  A.  The 
technique  allows  the  conversion  of  uncertainty  ellipses  to 
radii  of  circles  of  equivalent  probability.  Values  froa 
this  coaputation  for  the  radii  of  circles  of  90%  probability 
at  different  ranges  used  in  the  survey  area  for  15  hay  are 
displayed  in  figure  6.2  along  with  the  circles  of  90% 
probability,  the  aagnitude  of  the  radii  at  the  ranges  used 
for  the  study  do  not  exceed  5.4  aeters.  One  goal  of  this 
study  was  to  aeasure  the  difference  in  positional  accuracy 
between  AHTEMis  and  GPS.  A  aaxiaua  error  of  1.0  aeters,  an 
order  of  aagnitude  smaller  that  the  tested  system  accuracy, 
is  desirable.  At  the  present  time,  however,  there  are  few 
systems  that  can  provide  the  order  of  accuracy  system  for 
the  position  of  a  moving  vessel,  so  the  5.4  aeter  error 
observed  at  4000  aeters  had  to  suffice. 

Error  caused  by  reflection  of  the  signal  from  nearby 
objects  is  difficult  to  aeasure.  To  minimize  this  problem, 
the  offset  position  for  the  A8TENIS  unit  was  plaaned  so  as 
to  be  clear  of  possible  interference  froa  the  wooden 
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ARTEMIS  System  Error 


structure  and  the  equipaent  vehicle.  So  signal  interference 
was  observed  despite  the  ship*s  aotion.  This  source  of  data 
unreliability  was  therefore  not  considered  to  be 
significant. 

2.  ARTBHIS  System  Errors 

Errors  unique  to  the  ARTEMIS  systea  include  a 
pointing  error  source  and  the  lack  of  resolution  of  the 
aziauth  measurement  in  the  dynaaic  aode  of  operation. 

The  operation  of  the  ARTEMIS  systea  requires  a 


manual  alignaent  of  the  fixed 

antenna 

so 

that  it 

is 

perpendicular 

to  a  known  object. 

This 

is 

necessary 

to 

establish  an 

original  reference 

direction 

upon  which 

all 

angular  neasureaents  are  based.  It  is  possible  that  a 
pointing  error  could  be  produced  by  the  operator  if  the 
crosshairs  of  the  telescope  used  to  align  the  antenna  are 
not  exactly  sighted.  This  would  introduce  a  systematic 
error  into  all  subsequent  angle  aeasureaents.  To  ainiaize 
this  source  of  error*  the  alignaent  of  the  antenna  was 
checked  before  and  after  the  period  of  operation.  The 
original  aziauth  displayed  on  the  antenna  unit  display 
agreed  exactly  with  the  final  aziauth  value.  Any 
contribution  to  ship  position  error  was  considered  ainiaal. 


The  ARTEMIS  was  operated  using  the  dynaaic  node 
option  (see  Equipment  Descript ion--ARTEMIS) .  The  resolution 
of  the  systea  in  this  aode  is  1  aeter  for  displayed  range 
and  0.1°  for  aziauth  values.  This  resolution  Unit 

.  O 

introduces  an  uncertainty  of  _0.5  meters  in  range  and  0.05 
in  angle  measurements.  The  computed  uncertainty  ellipse  as 
expressed  by  radii  of  circles  of  equivalent  probability,  is 
shown  in  figure  6.2.  The  system  error  and  all  other  errors 
aentioned  above  are  assuaed  to  occur  within  the  scope  of 
this  resolution  uncertainty. 

Another  source  of  system  error  has  been  proposed  in 
a  study  conducted  in  June,  1980,  by  NOAA.  [Ref.  *9]  The 
test  compared  the  ARTEMIS  systea  with  two  ranges  measured  by 
two  Del  Norte  Trisponders  (a  short  range  microwave 
range/range  positioning  systea)  and  two  lines  of  direction 
measured  by  two  shore-based  theodolites.  The  results 
described  an  underway  bias  of  -0.05°  of  unknown  origin  in 
the  angle  aeasureaents.  This  discrepancy  has  not  been 

resolved.  The  preliminary  test  conducted  in  Monterey  in 
April,  1981  coapared  the  ARTEMIS  systea  with  three  lines  of 
direction  deterained  by  three  theodolites.  The  purpose  of 
the  test  was  to  resolve  the  underway  bias  problem.  Results 
are  enclosed  in  a  separate  report.  [Ref.  50] 
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8.  CONTROL  NETWORK 

The  entire  survey  was  oriented  relative  to  the  TRANSIT 
position  on  the  Coast  Guard  Pier  (GE0CEI7BB  ST  A  31370). 
Therefore,  errors  inherent  ia  the  position  of  this  station 
mark  affected  all  results.  Assuming  that  the  geographic 
coordinates  published  by  DR A/HTC  represented  the  best 
estimate  available  of  the  GEOCEIVEB  STA  31370,  the  recaining 
significant  sources  of  error  affecting  all  results  were 
survey  errors.  The  control  network  for  the  study  was 
determined  froa  the  traverse  and  it  was  assumed  that  all 
positions  used  in  the  study  were  oriented  relative  to  each 
other  on  the  BGS-72  ellipsoid.  A  change  in  the  ellipsoidal 
surface  within  the  area  encoapassed  by  the  traverse  was 
assumed  to  be  insignificant. 

Survey  errors  stea  froa  the  methods  used  for  obtaining 
the  initial  azimuth  for  orientation  of  the  survey,  the  type 
of  traverse  chosen,  and  the  geometry  of  the  traversed 
triangle. 

To  obtain  a  starting  azimuth  for  the  traverse,  two  sets 
of  observations  on  Polaris  were  recorded  at  the  SE0CEI7ER 
STA  31370.  A  striding  level  attachment  to  the  theodolite 
was  not  used.  The  instrument  level  bubble  was  observed 


repeatedly  daring  the  observations  to  ensure  constant 
instrusent  alignsent  with  the  vertical.  Failure  to  keep  the 
instrument  level  could  have  introduced  a  *  30"  uncertainty 

into  the  starting  azisuth  of  the  survey.  (Bef.  51]  To 
consider  the  significance  of  such  an  error,  all  stations 
were  recomputed  using  a  deviation  of  +  15"  from  the  computed 
azimuth.  The  resulting  shift  in  the  positions  of  MUSSEL  and 
USE  hoh  are  listed: 

OSS  SON  0.133  e  distance  at  230° 36*47.  16n  azimuth 

aosSEL  0.132  a  distance  at  230° 33*  4.08"  azimuth 

The  change  is  considered  insignificant  relative  to  the  other 
errors  described  above. 

The  survey  was  a  closed  loop  traverse  using  SEOCEIVEH 
STA  31370  as  the  origin  and  closing  position.  The  weakness 
in  this  method  involves  the  measured  distances.  A  method  of 
checking  the  observed  angles  in  the  triangle  exists  via  the 
angular  closure  if  a  significant  deviation  in  one  of  the 
angles  can  be  detected.  The  distance  measurement,  however, 
can  produce  an  undetected  bias  in  the  survey  if  the  distance 
measuring  instrument  contains  a  systematic  error.  This 
error  would  shift  all  the  surveyed  positions.  Correct 


measurement  of  distance  by  the  MBA  V  Telluroaeter, 
therefore,  was  important.  The  equipment  was  calibrated  over 
a  known  baseline  and  the  errors  were  within  the  centiseter 
range.  This  error  was  therefore  negligible. 

The  geoaetry  of  the  triangle  was  not  optiaua.  Ideally, 
the  triangle  would  have  been  equilateral  to  ainiaize  the 
source  of  observing  error.  The  choice  of  traversed  points 
was  based  on  the  location  of  preexisting  station  aarks  and 
the  desired  geoaetry  of  the  theodolite  fixes  to  the  ship. 
The  closures  computed  for  the  traverse  aet  third  order 
specifications  for  accuracy,  and  the  survey  was  not  adjusted 
to  reaove  the  closure  error.  [Ref.  52,  53J 

The  survey  errors  were  considered  to  be  within  the 
uncertainty  of  the  TRAHSIT  station  location  and  were 
therefore  not  considered  significant  for  the  study. 

C.  THEODOLITE  RET  BO  RK 

The  errors  generated  by  the  theodolite  network  are 
divided  into  equipment  errors,  operator  related  errors  and 
systea  errors,  of  which  the  latter  is  the  aost  significant. 

The  equipaent  used  for  the  theodolite  network  were  three 
1-second  theodolites  which  have  a  measurement  resolution  of 
i0.5w.  All  theodolites  were  in  good  condition. 


Instrunental  errors  of  a  theodolite  are  listed  is  Table  71. 
£ Hef .  54] 

The  aajor  source  of  equip aent  error  was  introduced  by 
not  changing  face,  the  process  of  observing  angles  a  second 
tine  with  the  telescope  reversed.  Due  to  the  nature  of  the 
test  coupled  with  1-ninute  fix  intervals,  it  was  not 
possible  to  observe  all  the  fixes  in  both  the  direct  and 
reversed  nodes.  As  an  exaaple  of  the  nagnitude  of  this 
error,  the  horizontal  colliaation  errors  were  deterained  for 
each  theodolite  prior  to  the  study.  The  results  are  listed 
below: 

Xhsaflflllto  Sail  iaaU<?n  sclqj 

*  14405  -  1.0" 

*  14482  ♦  2.7" 

*  14452  -  0.7" 

All  the  aeasured  angles  were  adjusted  for  this  error  as 
described  in  the  processing  section.  The  largest  influence 
on  the  reaaining  instruaental  errors  was  the  non- vert icality 
of  the  aain  axis.  If  the  non-verticality  were  large  (for 
exaaple,  one  division  off-level  in  the  lower  spirit  level 
bubble) ,  the  angle  aeasured  could  be  in  error  by  7". 
[8ef.  55]  To  ainiaize  this  error,  each  theodolite  was 
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TABLE  VI.  Theodolite  Instrument  Errors 


ERROR 

Eccentricity,  misplaced  indices 

Graduation  of  circles 
Micrometer  run  and  vernier  error 

Non-verticality  of  main  axis 

Bad  setting  of  the  upper  spirit 
level 

Vertical ity  of  cross  wires 

Horizontal  collimation 

Vertical  collimation  and  index 
setting  error 


ELIMINATED  BY: 


Reading  both  sides  of  circle 
and  changing  face 

Changing  zeros 

Selection  of  suitable  zero 
setting 

Error  is  not  cancti  '■’$#  4ny 
observi ng  procedure 

Correct  setting  nr  na¬ 

tion  of  correction 

Changing  fac*  same 

part  of  wire 

Changing  face 

Changing  face 
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relabelled  if  necessary  between  sets  of  observations.  The 
aagnitude  of  this  and  the  colliaation  errors  were  considered 
to  be  insignificant  relative  to  the  systea  errors  discussed 
later  in  this  section. 

The  operator-related  errors  are  a  function  of  the 
individual  perforaance  of  the  observers.  The  experience  and 
personal  coafort  of  the  observer,  for  example,  are  two  aajor 
factors  which  can  influence  the  observed  angles.  An 
observer's  perforaance  and  resultant  iapact  on  the  quality 
of  the  results  depend  on  the  judgaent  and  reaction  tiae  of 
the  observer  and  are  influenced  by  variables  such  as  fatigue 
and  age.  The  subjective  nature  of  the  problea  makes  an 
estimate  of  the  error  difficult,  therefore,  for  this  study 
the  error  was  considered  random. 

The  aajor  source  of  error  in  the  theodolite  positioning 
systea  dealt  with  the  tiaing  of  the  fix.  A  delay  existed 
between  the  tiae  the  'aark'  was  called  over  the  radio  and 
the  actual  tiae  the  angle  was  observed.  The  tiaing  problea 
was  most  critical  when  the  ship  was  close  to  one  of  the 
three  observer  locations  because  the  rate  of  change  of  the 
bearing  of  the  ship  was  greatest  at  this  tiae.  Hence  a 
saall  delay  generated  a  large  aaount  of  error  in  the  angle 
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recorded.  To  indicate  the  effect  this  error  source  say  have 
had  on  the  position  of  the  ship,  tvo  fixes  vere  chosen 
arbitrarily — one  at  close  range  to  each  of  MUSSEL  and 
6S0CEITEB  ST1  31370.  A  one  ninute  value  was  added  to  the 
angle  which  vas  aost  affected  by  the  timing  error  and  the 
fixes  vere  recomputed.  The  same  process  vas  repeated  for 
the  three  fixes  with  the  one  minute  value  subtracted  from 
the  critical  angle  in  each  fix.  The  results  follov: 

Humber  Original  Fix  *  1  Minute  -  1  Minute 

210  36°37*20. 62282"  36° 37* 20.6 146 8"  36°37* 20.63081" 

121°53*41. 97184"  121°  S3* 4  1 . 9693 3"  1 2 1° 53 •  4 1 . 97426" 
259  36°36'45. 88048"  36°  36*45.87574"  36 °36 • 45. 89082« 

121°52*54. 05062"  121  °  52* 54. 04635"  121°52* 54. 05459" 
Ho  significant  error  vas  generated  from  this  source. 

Another  source  of  system  error  vas  associated  vith  the 
geometry  of  the  fix.  Ideally,  the  three  theodolite  angles 
would  describe  tvo  equilateral  triangles  to  the  ship; 
however,  within  the  operating  area,  the  geometry  varied.  As 
an  indication  of  the  size  of  this  error,  representative 
fixes  vere  chosen  and  the  circular  map  accuracy  standards 
(CHAS)  of  each  of  the  combinations  was  determined.  Figure 
6.3  illustrates  the  approximate  error  caused  by  variations 
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Fig.  6.3.  Representative  Theodolite  Intersections  for  Error  Ellipse  Computations 


in  the  geometry.  iithin  the  area  outlined  in  the  figure  the 
geometry  creates  a  CHAS  of  not  lore  than  10.4  aeters.  This 
is  therefore  the  source  of  the  largest  probable  error  in  the 
theodolite  oethod  of  positioning.  Nevertheless,  this  method 
in  cosbination  with  the  ARTEMIS  system  provided  the  best 
estimate  of  the  ship  location. 

0.  TRANSIT  POSITIONING  SYSTEM 

The  aethod  used  to  deteraine  the  TRANSIT  position 
(GEOCEIVBB  STA  31370)  is  the  aost  accurate  satellite 
positioning  technique  presently  available.  Error  sources 
are  reduced  both  during  data  collection  and  during 
postprocessing  of  the  recorded  values. 

JMA  places  requireaents  on  the  collection  of  data  that 
will  ainiaiza  error.  For  this  study,  only  satellite  passes 
whose  aaziaua  altitude  ranged  between  10°  and  30°  were 
recorded.  Meteorological  data  (teaperature,  pressure,  and 
humidity)  were  collected  for  each  pas3,  and  a  determination 
of  the  receiver  clock  drift  was  computed  after  each  pass. 
The  amount  of  clock  drift  as  a  function  of  time  gave  a 
real-time  indication  of  receiver  performance.  Equipment 
checks  were  conducted  routinely,  and  the  results  were 
included  with  the  data  package  submitted  to  DMA . 
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An  error  budget  of  the  TRANSIT  systea  for  point 
positioning  is  listed  in  Table  VII.  [Ref.  56]  The  precise 
epheaeris  aethod  of  postprocessing  accounts  for  all  of  these 
errors:  1)  for  each  pass  the  tropospheric  correction  is 
coaputed  froa  the  aeteorological  data;  2)  the  ionospheric 
correction  is  deterained  by  observing  the  shift  in  frequency 
between  the  two  transaitted  frequencies;  3)  the  receiver 
clock  error  is  known;  4)  the  satellite  orbits  are  observed 
by  the  tracking  stations  so  that  precise  satellite  positions 
are  known;  5)  errors  due  to  uncertainties  in  the 
geopotential  aodel  and  surface  drag  forces  are  accounted  for 
during  coaputation;  and  6)  the  receiver  altitude  relative  to 
the  ellipsoidal  surface  is  coaputed  approximately  froa  the 
given  MSL  elevation. 

The  accuracy  of  the  coaputed  TRANSIT  position  is  stated 
as  1.5  aeters  per  axis  within  a  90%  confidence  interval. 
This  is  assuaed  to  provide  the  best  estiaate  possible  of  the 
GEOCEIVER  STA  31370  relative  to  the  WGS-72  coordinate 
systea. 

S.  GPS  ERROR  SOURCES 

Global  Positioning  Systea  error  sources  are  separated 
into  two  aain  categories,  the  space  vehicle  systea  errors 


TABLE  VI I. TRANSIT  System  Errors 


SOURCE  ERROR  ( 


Uncorrected  propagation  effects 
(ionospheric  and  tropospheric  effects)  1-5 

Instrumentation  and  measurement  noise 
(local  and  satellite  oscillator  phase 
jitter,  navigator's  clock  error)  3-6 

Uncertainties  in  the  geopotential  model 
used  in  generating  the  orbit  10-20 

Incorrectly  modelled  surface  forces 
(drag  and  radiation  pressure  acting  on 
the  satellites  during  the  extrapolation 
interval)  10-25 

Uncertainties  in  navigator's  altitude 
(results  in  a  bias  in  longitude)  10 

Ephemeris  rounding  error  5 
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and  the  GPS  user  receiver  errors.  There  are  seven  possible 
sources  of  error  in  the  space  vehicle  systea — errors 
associated  with  satellite  equipeent,  satellite  clock, 
epheaeris,  ataospheric  delays,  signal  aultipath  effect, 
satellite  update,  and  satellite  propulsion  firing.  Table 
III  lists  these  error  sources  and  the  approxiaate  user 
equivalent  range  error  for  each.  [fief.  57] 

The  satellite  equipaent  error,  also  teraed  the  group 
delay,  is  caused  by  differences  and  uncertainties  in 
processing  and  transaission  of  the  aessage  through  the 
individual  satellite's  circuitry.  Bach  satellite  has  the 
group  delay  calibrated  prior  to  launch  and  this  known  group 
delay  is  included  in  the  satellite  tiae  corrector  within  the 
navigation  aessage.  However,  unknown  delays  in  signal 
processing  and  circuitry  not  accounted  for  by  the  corrector 
contributes  an  estiaated  one  aeter  to  the  error  budget. 

Satellite  clock  error  is  caused  by  variations  in  each 
satellite  clock.  Satellite  clock  tine  may  depart  by  up  to 
976  microseconds  froa  the  correct  GPS  tiae.  Corrections  for 
the  predicted  clock  variations  are  included  in  the 
navigation  aessage.  Residual  clock  variations  contribute  to 
the  error  budget  and  are  included  in  epheaeris  errors  since 


the  clock  errors  are  saall  and  have  the  saae  effect  as 


epheaeris  errors. 

Epheaeris  errors  are  the  result  of  satellite  clock 
variations,  drift  by  the  aonitor  station  clocks,  aonitor 
station  signal  processing  delays,  lack  of  precise  aonitor 
station  geodetic  positions,  solar  wind  and  pressure 
paraaeter  variations,  and  earth  gravity  model  error.  These 
sources  all  contribute  to  the  epheaeris  error  budget.  The 
epheaeris  error  budget  is  aonitor ed  over  long  tiae  periods 
and  translated  by  the  Master  Control  Station  into  an 
User-Equivalent  Range  Error  (ORE),  when  the  ORE  exceeds  the 
prescribed  value  of  4  aeters  the  satellite  epheaeris  data  is 
uploaded  to  reduce  the  ORE.  (Ref.  58]  Most  of  the  ORE  is 
included  into  the  user  clock  corrector  obtained  in  the  GPS 
position  coaputation.  However,  not  all  of  the  epheaeris 
error  is  reaoved  during  the  position  coaputation,  and  the 
residual  error  has  bees  estiaated  to  be  approxiaately  1.5 
aeters.  See  Table  III. 

ktaospheric  delays  in  signal  propagation  are  the  result 
of  two  effects,  ionospheric  delay  and  tropospheric  delay. 
Of  these,  the  ionospheric  delay,  which  depends  on  the 
frequency  of  the  transaitted  signal,  is  greater.  With  the 
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two  frequencies,  LI  and  L2,  transmitted  by  each  satellite  to 
the  user  receiver,  the  ionospheric  propagation  delay  can  be 
modelled.  rising  this  model,  the  majority  of  the  propagation 
delay  is  accounted  for,  but  a  residual  delay  results  in  an 
error.  The  other  part  of  the  atmospheric  delay  error 
depends  on  the  tropospheric  propagation  delay.  Propagation 
delay  in  the  troposphere  is  a  function  of  signal  path  length 
through  the  layer,  the  receiver  air  temperature,  and  air 
pressure  and  humidity.  Since  the  tropospheric  delay  is  not 
a  function  of  frequency  and  is  highly  localized,  the 
tropospheric  delay  corrector  has  not  been  incorporated  into 
the  atmospheric  correction  computation  and  is  therefore  an 
additional  source  of  error.  These  two  error  sources  combine 
to  yield  an  atmospheric  delay  error  estimated  to  have  a 
magnitude  of  between  2.4  to  5.2  meters. 

Ionospheric  scintillation  of  the  GPS  signal,  in  addition 
to  the  ionospheric  propagation  delay,  cause  errors  in 
positioning.  Scintillation  has  the  effect  of  random 
interference  with  the  propagated  signal  as  the  signal  is 
transmitted  through  the  upper  ionosphere.  [Ref.  59] 
Interference  caused  by  scintillation  combines  with  the  GPS 
user  receiver  processing  noise  to  cause  a  random  ranging 
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error  and#  if  the  interference  is  severe#  interruption  of 
the  pseudo-random  number  binary  code  causes  loss  of  signal. 
Errors  and  signal  losses  caused  by  scintillation  cannot  be 
predicted  or  estiaated. 

Signal  aultipath  errors  result  froa  the  reception  of  a 
satellite's  signal  transnitted  via  different  propagation 
paths.  Received  signals  thus  say  have  been  distorted  or 
reduced  in  strength  to  a  sufficient  degree  so  as  to  affect 
the  range  determination  from  the  satellite.  Error  due  to 
aultipathing  is  highly  localized  and  cannot  be  modelled. 
Estimates  of  the  error  magnitude  are  between  1.2  and  2.7 
meters. 

The  last  two  space  vehicle  system  error  sources  are 
attributable  to  the  flaster  Control  Station  (MCS)  satellite 
accuracy  and  maintenance  functions.  Usually  the  NCS  uploads 
each  satellite  at  a  given  time  each  day  and  the  ephemeris 
data  is  changed(Table  7III).  should  the  user  not  be  aware  of 
occurrence  of  upload  and  if  the  receiver  therefore  does  not 
obtain  the  new  ephemeris  data#  ranging  to  the  satellite  in 
question  will  contain  a  possibly  very  large  error.  In 
addition#  if  the  HCS#  through  the  Beaote  Stations# 
determines  that  the  four  skewed  reaction  wheels  cannot 
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maintain  satellite  antenna  alignment  with  earth, 
hydrazine  propulsion  systea  is  activated  and  the  satellite 
is  rotated.  Once  the  propulsion  system  has  used,  the 
satellite  orbit  cannot  stabilize  immediately,  and  the 

epheneris  data  will  be  inaccurate.  k  period  of  up  to 
several  hours  may  be  required  for  orbit  stabilization  and 
for  the  new  epheaeris  data  to  be  computed  and  uploaded. 
During  that  period,  the  ranging  computations  to  the 
satellite  will  have  large  errors.  The  magnitude  of  these 
two  error  sources  cannot  be  estimated  but  is  assumed  to  be 
several  orders  of  magnitude  larger  than  the  previously 
stated  errors. 

The  second  error  category  is  that  related  to  the  GPS 
user  receivers.  These  sources  of  error  associated  with  the 
user  GPS  receiver,  the  Manpack,  are  receiver  signal  noise, 
receiver  vessel  velocity,  signal  resolution,  method  of 
signal  processing,  epheaeris  update  capability,  and,  for  the 
differential  GPS  computations,  system  accuracy  comparison. 

Deceived  signal  noise  errors  result  from  interference 
introduced  into  the  receiver  GPS  signals  through  processing 
and  internal  receiver  circuitry.  Additionally,  quality  of 
signal  resolution  is  dependent  upon  receiver  type  and  can  be 


a  source  of  error.  The  coebined  errors  froa  these  two 
sources  are  estimated  to  be  1.S  aeters.  [Ref.  60] 
Estimates  of  these  error  sources  for  the  8V0E  have  not  been 
aade  but  are  estinated  to  be  of  the  saae  magnitude. 

Telocity  of  the  receiver  vessel  results  in  GPS  errors  of 
varying  magnitude.  The  higher  the  velocity,  the  larger  the 
error.  This  error  is  due  to  the  time  lag  required  to  obtain 
and  process  each  N AVSTAR  satellite  range  with  the  mvoe 
single  channel  receiver,  to  compute  the  average  vessel 
velocity  and  to  use  the  average  velocity  for  positional 
computations.  The  (Unpack  requires  a  six  second  time 
interval  to  obtain  and  process  the  pseudo-range  from  each 
satellite  used.  Thus  the  fix  position  is  an  extrapolation 
of  the  four  satellite  ranges  to  obtain  a  convergent  fix 
position.  At  one  knot,  the  vessel  travels  0.51  aeters  in 
one  second.  During  the  full  twenty-four  seconds  required  to 
obtain  four  satellite  ranges,  the  vessel  would  have 
travelled  12.24  aeters.  With  the  vessel  travelling  at  four 
kncts,  the  distance  travelled  during  the  total  time  interval 
required  for  a  GPS  fix  position  would  be  48.96  meters. 
Processing  the  pseudo-ranges  froa  four  satellites  requires 
the  solution  including  the  average  vessel  velocity  during 
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the  24-second  ranging  period.  Use  of  this  average  neglects 
the  relative  velocity  of  the  receiver  antenna  which  is 
■ousted  on  the  last.  During  this  test,  the  GPS  antenna  was 
aounted  9.06  aeters  above  the  waterline  of  the  R/7  ACANIA. 
Assuaing  five  degrees  of  pitch  and  roll  and  a  vessel  roll 
period  of  IS  seconds,  relative  average  antenna  velocity  was 
,0.16  aeters  per  second.  The  antenna  velocity  error 
translates  into  a  random  error  of  0.16  aeters.  The  other 
velocity  errors  have  not  been  estimated  but  are  assuaed  to 
be  less  than  one  neter  for  the  low  vessel  velocities  used  in 
the  test. 

one  error  source  mentioned  previously  requires  further 
clarification:  error  due  to  the  relative  position  at  the 

tiae  of  the  fix.  The  Nanpack  requires  up  to  twenty-four 
seconds  to  obtain  and  process  the  satellite  ranges.  The 
position  of  the  fix  and  the  associated  tiae  of  fix  are  then 
recorded.  However,  the  tiae  of  the  fix  may  not  have  been 
the  saae  time  the  computational  results  were  printed.  Thus 
the  output  position  may  not  have  been  correlated  with  the 
correct  theodolite  computed  ship  position  since  the 
theodolite  lines  of  position  were  observed  when  the  GPS 
position  was  printed.  Errors  froa  this  bias  source  may 


range  froa  12.24  eaters  at  one  knot  of  vessel  velocity  to 
48.96  aeters  at  four  knots.  This  error  varied  with  the 
processing  time  required  to  obtain  a  convergent  position, 
in  estiaate  of  this  error  aagnitude  was  not  aade  because  of 
its  vide  variability. 

The  type  of  receiver  and  the  aethod  of  signal  processing 
also  contributed  to  the  error.  Bultichannel  receivers 
designed  for  Phase  I  testing  have  a  lover  aagnitude  error 
for  satellite  signal  reception  and  processing  than  do  the 
single  channel  receivers,  such  as  the  Hanpack.  In  the 
aultichannel  receivers,  four  or  aore  satellite  pseudo-range 
data  are  taken  siaultaneously .  The  single  channel  receiver 
sequentially  receives  each  satellite  and  coaputes  the 
pseudo-range  data.  Siaultaneous  data  reception  reduces  the 
user  velocity  error  and  therefore  has  less  position  aliasing 
than  that  found  in  the  single  channel  receivers. 

The  aethod  of  3ignal  processing  aay  induce  positional 
errors  through  the  use  of  non-optiaal  weighting  functions  in 
the  position  computation.  Should  the  epheaeris  or  ranging 
data  change  significantly,  as  happens  during  uploading, 
signal  processing  using  Kalaan  filters  as  used  in  the 
danpack  do  not  change  the  coaputed  position  rapidly  enough 


to  reflect  the  large  changes  in  ranges  and  result  in  a  slow 
change  drift  of  the  GPS  position.  Estimated  tine 
reguireaents  for  user  receivers  to  stabilize  after  an  upload 
are  up  to  two  hours.  [Ref.  t  i  ] 

The  capability  of  the  user  GPS  receiver  to  receive 
epheneris  data  iaaediately  after  upload  reduces  the 
epheaeris  errors  caused  by  using  an  outdated  epheoeris.  The 
epheaeris  data  is  repeated  every  thirty  seconds  and  its 
length  reguires  sections  of  two  message  subframes. 
[  Ref .  62  .  63Q.  The  use  of  current  epheaeris  data  in 
the  computations  reduces  the  error  for  the  epheaeris  to  the 
value  shown  in  table  III.  The  Ran  pack  has  the  capability  to 
accept  new  epheaeris  data  iaaediately,  but  during  the  test 
large  errors  were  introduced  into  the  coaputation  on 
satellite  epheaeris  upload  because  this  selected  optional 
function  either  did  not  operate  properly  or  was  neglected  by 
the  operators. 

For  the  differential  GPS  user,  an  additional  error 
source  is  introduced  into  the  error  budget:  the  assumption 
that  both  GPS  receivers  will  observe,  obtain  the 
pseudo- ranges ,  and  compute  the  position  solution  in  exactly 
the  saae  Banner.  If  this  were  true  then  both  receivers 
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should  be  able  to  coopute  the  identical  position  if  placed 
at  the  saae  position  at  the  saae  tiae.  With  the  Manpacks, 
the  assuaption  that  the  units  were  identical  and  the 
position  solutions  were  identical  was  not  aade  because  prior 
testing  established  that  the  receivers  differed  in  position 
coaputations  by  as  auch  as  3  08  aeters.  [Ref.  64] 

Por  the  static  testing  on  11  May,  one  error  source  not 
previously  discussed  was  the  vertical  separation  of  the  two 
Manpack  antennas,  colocated  over  the  GEOCEIVER  STA  31370 
control  station,  and  the  possible  interference  of  the 
upperaost  antenna,  Systea  1,  with  the  GPS  signals  received 
by  the  lower  Systea  2  antenna.  Separation  of  the  two 
antennas  was  2.1  aeters.  Calculations  using  the  speed  of 
propagation  in  a  vacuua  and  the  two  GPS  frequencies  produce 
wave  lengths  of  0.19  aeters  and  0.24  aeters.  Since 
separation  of  the  antennas  was  auch  greater  than  either 
wavelength  theoretically  there  should  have  been  no 
positional  errors  due  to  antenna  interaction  or 
interference.  Differing  elevations  of  the  antennas  above 
the  ellipsoidal  surface  also  aay  have  induced  errors  in  the 
position  coaputation.  Sstiaated  errors  for  the  elevation 
differences  are  less  than  0.1  aeters.  Total  errors  for  the 
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antenna  vertical  elevation  separation  and  antenna 
interference  or  interaction  have  not  been  estiaated. 

Thus,  the  error  budget  for  the  Global  Positioning  System 
suggests  a  combined  error  of  up  to  57.35  meters.  Using  the 
differential  technique,  the  bias  errors  can  be  reduced  in 
magnitude  or  eliminated.  These  bias  errors  include 
ionospheric  propagation  delays,  tropospheric  propagation 
delays,  variations  from  predicted  epheaerides,  and  satellite 
clock  perturbations.  Teasley,  Hoover,  and  Johnson  (1980), 
evaluated  the  error  budget  for  a  more  advanced  GPS  receiver 
and  found  the  predicted  filtered  solution  error  was  reduced 
from  approximately  11  meters  to  approximately  3  meters  with 
differential  GPS  processing.  [Ref.  65]  Error  budgets  for 
the  Hanpack  and  the  differential  correction  vector 
processing  method  used  in  this  test  did  not  attain  the  four 
meter  accuracy  but  were  not  expected  to  because  the 
receivers  used  were  less  accurate,  and  of  different  types. 


134 


?II.  RESOLTS 


Results  of  the  GPS  tests  have  been  separated  into  four 
sections:  stationary  GPS  evaluation,  evaluation  of  non- 
differential  GPS,  differential  GPS  evaluation  with  static 
correction,  and  evaluation  of  differential  GPS  without 
static  correction.  Tables  IV  and  VIII  list  the  optimum 
satellite  observation  periods  and  the  N AVSTAR  upload  times 
for  each  day  of  the  testing  period. 

A.  STATION  ART  GPS  EVALUATION 

Stationary  GPS  receiver  evaluation  was  separated  into 
the  stability  tests,  the  static  tests  and  the  tests  of  the 
stationary  GPS  receiver  used  during  differential  testing. 
Pigures  7.1,  7.2,  and  7.3  illustrate  the  results  for  the 
stability  test  conducted  on  9  and  10  Hay.  The  receiver  was 
less  stable  than  expected  and  depended  to  a  great  extent  on 
the  MAVSTAR  uploading  times  and  the  ORE  of  NAVSTAR  1  (fig. 
7.4  and  7.5) .  Correlation  between  times  for  the  NAVSTAR 
upload  and  changes  in  the  stability  can  be  observed  at  0620 
on  9  Hay  (fig.  7.1),  0655  on  9  Hay  (fig.  7.2),  0815  on  10 
Hay  (fig.  7.3),  and  0845  on  10  Hay  (fig.  7.3).  The 
relationship  between  stability  and  ORE  is  illustrated  by 
comparison  of  the  ORE  for  NAVSTAR  1  at  0640  on  9  Hay  (fig. 
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LONGITUDE  _  LATITUDE 


Figure  7.4. 


GREENWICH  MEAN  TIME  (HOURS) 


7.4)  with  the  stability  plot  on  figure  7.1  for  the  sane  tine 
period.  On  figure  7.2  the  stability  plot  variation  at  0745 
on  10  Hay  renains  unexplained  but  nay  have  been  due  to  a 
receiver  naif unction. 

Sesults  for  tests  on  11  Say  (fig.  7.6,  7.7,  and  7.8), 
exhibited  less  stability  than  expected.  Also,  systen  2  vas 
less  stable  than  systen  1.  The  largest  difference  in 
reading  between  the  two  receivers  occurred  at  0750  (fig. 
7.7).  This  difference  nay  have  been  caused  by  the  location 
of  the  antenna  of  systen  2  under  the  antenna  of  systen  1. 
At  approxinately  0800,  the  satellite  constellation  elevation 
angles  were  the  largest,  between  56  °  to  80°  above  the 
horizon  (fig.  7.9) ,  and  antenna  interference  nay  have  led  to 
the  large  differences.  Neither  the  satellite  constellation 
azinuth  angles  (fig.  7.10)  nor  the  user  range  errors  (fig* 
7.11)  appeared  to  have  any  effect  on  the  stability. 

Conparison  of  the  stationary  OPS  receiver  positions  with 
the  position  obtained  with  the  geoceiver  is  depicted  in 
figure  7.12  for  the  stability  and  static  tests.  The 


conparison  indicates  that  there  is  a  nearly  constant  offset 
of  the  two  positions  with  the  GPS  position  almost  due  south 
of  the  geoceiver  position.  Further  comparisons  between  the 
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Fig.  7.3. 
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GREENWICH  MEAN  TIME  (HOURS) 
Fig.  7.9. 
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Figure  7.11 . 
TOTAL  USER  RANGE  ERROR 


GREENWICH  MEAN  TIME  (HOURS) 


SEQCEIVER  VS.  GPS  POSITIOM  COMPUTATION 

GEOCEIVER  POSITION  VS.  GPS  SHORE-BASED  POSITION 
FROM  OATA  COLLECTED 
9,  10,  II  MAY  STATIC  AND 
STABILITY  TESTING 
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stationary  GPS  receiver  position  and  the  geoceiver  position 
reinforce  this  conclusion.  In  figures  7.13,  7.14,  and  7.15, 
the  aean  distances  between  the  two  positions  are  127,  173 
and  143  aeters,  respectively.  The  higher  aeans  for  the  data 
shown  in  figures  7.14  and  7.15  are  due  to  the  UBS  (fig.  7.16 
and  7.17)  for  those  dates  and  the  NAVSTAR  upload  frequency. 
If  the  data  froa  13  Nay  is  deleted  because  of  the  high  DBS 
for  that  date  (fig.  7.16),  and  the  aedian  distance  used  for 
15  Nay  is  taken  to  be  a  representative  value  (130  aeters  on 
fig.  7.15),  The  difference  is  approxiaately  123  aeters 
Between  the  two  positions.  The  aziauth  was  also  nearly 
constant  with  the  observed  median  values  of  18 f  ,  207°  and 
178 3  froa  south  as  shown  in  figures  7.18,  7.19,  and  7.20 
respectively.  Neglecting  the  data  on  13  May  because  of  ORE 
values  (fig.  7.16),  the  reaaining  two  values  are  nearly 

O 

identical  and  have  a  aean  of  179.5  froa  south.  Thus  the 
data  indicates  a  nearly  constant  offset  between  the  position 
obtained  with  the  geoceiver  and  the  position  obtained  with 
GPS.  This  difference  is  about  130  aeters  and  has  an  almost 
north-south  orientation. 
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Figure  7.14 
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FRECUENC IES  SAHRL6  SIZE 


TOTAL  USER  RANG 


GREENWICH  MEAN  TIME  (HOURS) 


FREQUENCIES  sahrle  sue 


mCUENCUS  SA1PLE  SUE 


Figure  7.19 


-! 


I 

♦ 

on 

i 

i 

♦ 

on 

i 

i 

♦ 

*n 

i 

i 


*  U»*  *  «#•*•»»•••»««  N 
»  U.«  %•#*-••*■**•#*•* 


u. 

**. 

u. 


I 

« 

♦ 
O  I 

I 


t«*rtttWX«t*«»’**ii«rt*r*r**r»»rr**rTrTTrttrTri 


U.  «•#«••*••••••••••••••••••••••♦••••••••••  ••••  *  •  ••• 


♦ 


UJUMMMMIUl 

sOir»cxrw\r»»- 

>0 

rO(>rs«C 

.  ■■ 


U*«UA 

sa? 


3  2222*^ 

«X«««2 

—3393333 


•••«••  ««•  ■> 


vt20indMNo« 

cx  •••••* 


_ 

kOOOOOO 
?♦♦ • I ♦♦ 

mimummiym 

z 

CNtPOWOfl 

INtTNOO0 

jfUNOMftNN 

<O<O'0W«O 

«•••••• 

2  I 


— O  «M  20^0* 

•-*-  X  :*►-■*< 

o  ujoc»*»* 

►—  — 

«/»»** 

24 

u*» 

HIM 


UAUA 

<«» 


30 

QUJ 

VlO 


oooooo 

♦  ♦  i  ♦  ♦  ♦ 
uMIUWlUUi 
W»P  Nrt 
0«* 

— **v*A«rv\>a 

0404«rt 

•  ••••» 


!2  oi25»  S 

Ul^  440  UJ&. 
ac—  uj— 
o*/»  «o*oua<zO 
uiO  4»-Oui4** 

30. 


NNNNNNN 


T— 


30 

XUJ  IUUMA«UUIUJEM 

MS 

«T  -.  OA^TO'fl— 

V-TK^CfO  —  EO— 
U^f«NOVH»4 


M.  I  -»  9 


OOjlN04>|fvO##M 

2 


4  ZZV9«utu 
«  2«Ztt 

£  <«UIUJ«V 

w?-VMtf 

440-ciog* 

IXIUJlAJCt— Ui«C 

uttHCXOX 


156 


FPECUENCIiS  SAMPLE  SUE 


Figure  7.20 
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3.  EVALUATION  OF  NONDIFFESSNTIAL  GPS 

Evaluation  of  the  non-differential  GPS  positions  was 
done  on  a  set  by  set  basis.  Pigures  7.21  through  7.31  depict 
the  distance  in  eeters  between  the  theodolite-computed  ship 
position  and  the  shipboard  GPS  receiver  position.  Sets  1,  2 
and  3  have  the  smallest  distance  difference,  as  shown  in 
figures  7.21,  7.22,  and  7.23,  with  mean  values  between  118 
and  133  meters  and  standard  deviations  ranging  from  7.45  to 
16  meters.  These  smaller  mean  distances  are  a  result  of  the 
low  user  range  errors  for  12  Nay  as  shown  on  figure  7.32. 
The  highest  mean  distances  occurred  during  sets  4,  9,  10, 
12,  and  13  when  the  UBE  values  were  high  and  the  NAVSTAB 
uploads  occurred  (fig.  7.16  and  7.17).  The  remaining 
distances  vary  between  144  and  278  meters. 

The  results  during  the  best  satellite  observation 
periods  are  depicted  in  figures  7.33,  7.34,  and  7.35,  where 
the  non-differential  GPS  ship  position  was  compared  to  the 
theodolite- computed  ship  position,  the  distance  between  the 
two  being  plotted  as  a  histogram.  Mean  distances  on  figures 
7.33  and  7.35  are  representative  of  the  five  satellites 
observed  and  the  improved  accuracy  of  the  system,  with  mean 
distances  of  134  and  149  meters  and  standard  deviations  of 
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FAECUEMCIES  SAMPLE  SUE 


Figure  7.21 
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MECUENC1ES  SAMPLE  SUE 


Figure  7  22 
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EAECUENCliS  SAMPLE 


Figure  7.23 
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fR£3U£NCUS  SAMPLE  SIZE 


Figure  7.24 
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FtlECUENCIES  SAMPLE  SUE 


Figure  7.25. 
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mcutuciES  sample  sue 


Figure  7.26 
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flltllUtfcCUS  sample  sue 


Figure  7.27 


fPECUEKCIES  SAMPLE  Si  EE 


mcufNCics  sample  sue 


Figure  7.31 
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Figure  7.33 
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10  and  52  aeters.  The  large  standard  deviation  for  15  May 
is  apparent  in  the  biaodal  eapirical  frequency  curve  (fig. 
7.35)  and  is  caused  by  the  SAT  STAS  upload  (fig.  7.17).  in 
this  case,  the  GPS  receiver  did  not  obtain  the  new  epheneris 
until  halfway  through  the  tiae  period  of  the  observations. 

4 

The  data  for  13  Hay  (fig.  7.34)  reflects  frequent  UAVSTAR  1 
uploads  and  the  resulting  high  USE  values  (fig.  7.16)  for 
that  date. 

C.  DIFFERENTIAL  GPS  EVALUATION  5JITH  STATIC  CORRECTION 

The  results  of  the  evaluation  of  the  data  of  12  May 
(sets  1,  2,  and  3),  computed  using  the  differential  GPS 
positions  and  the  static  corrections,  are  shown  in  figures 
7.36,  7.37,  and  7.38.  Mean  difference  between  positions  for 
the  theodolite  ship  position  and  the  GPS  position  are  29.1, 
27.4,  and  64.8  meters,  with  standard  deviations  of  7.7,  7.5, 
and  22.7  aeters.  These  distances  are  larger  than  those 
deterained  for  the  theodolite  ship  position  and  the 
differential  GPS  position  (fig.  7.39,  7.40  and  7.41).  This 
higher  aean  distance  indicates  that  the  static  correction, 
as  applied,  results  in  a  larger  error  in  the  differential 
GPS  position.  Therefore,  further  evaluation  of  differential 
GPS  with  static  corrections  were  not  conducted. 
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Figure  7.37. 
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Figure  ?.41 . 
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D.  EVALUATION  OP  DIFFERENTIAL  GPS  NITHOOT  STATIC 
COBBSCTIONS 

Evaluation  of  differential  GPS  without  the  static 
corrections  were  computed  by  sets  for  the  five  satellite 
constellations  for  each  day  of  testing.  Set  1  (fig.  7.39) 
shows  a  mean  distance  of  10.1  meters  between  the  two 
positions  and  a  standard  deviation  of  5.0  meters.  The 
median  azimuth  between  the  two  was  257°  from  south  (fig. 
7.42) ,  which  also  was  the  approximate  course  of  the  vessel 
during  the  test.  These  values  indicate  that  the  error  is 
caused  by  the  bias,  or  alias,  in  GPS  position  due  to  the 
Nanpack  24-second  fixing  interval  and  positioning  algorithm. 
For  this  set,  the  GPS  position  was  obtained  after  the 
theodolite  position. 

During  set  2  the  mean  distance  difference  between  the 
two  positions  was  11.7  meters,  and  the  standard  deviation 
was  6.8  meters  (fig.  7.40)  ;  the  median  azimuth  was  71° 
(fig.  7.43) .  These  values  were  slightly  greater  than  those 
of  set  1  due  to  the  increased  UBE  for  NAVSTAR  1  during  the 
set  (see  0700  on  fig.  7.32)  Additionally,  the  azimuth  was 
nearly  180°  from  that  computed  in  set  1,  with  the  vessel  on 


the  3ame  track  as  set  1.  Onlike  the  error  obtained  from  the 
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Figure  7.42 
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Manpack  bias  on  set  1 ,  the  error  in  set  2  was  due  to  a 
siailar  time  delay  bias  because  the  theodolite  position  was 
obtained  after  the  GPS  position. 

Set  3  values  on  figures  7.41  and  7.44  are  siailar  to 
those  of  set  2,  except  the  mean  distance  values  were 
greater — 44.7  meters  for  set  3  versus  11.7  meters  for  set  2. 
Mean  distance  increases  for  set  3  were  due  to  the  upload  of 
MAVSTAR  1  (fig.  7.32)  and  the  availability  of  only  a  four 
satellite  constellation.  Set  2  had  a  five  satellite 
constellation  available. 

?or  set  4 ,  figures  7.45  and  7.46  show  the  distances  and 
aziauths  for  the  set.  These  histograms  fail  to  Illustrate 
the  constantly  increasing  distances  between  the  two 
positions.  These  increasing  distances  result  froa  the 
uploads  of  NAVSTA8*  s  5  and  6  and  the  large  U RE  for  SAVSTAR  1 
(fig.  7.16).  Onboard  the  vessel  the  new  epheaeris  data 
uploaded  for  NAVSTAR  *s  5  and  6  was  not  obtained  until  late 
in  the  set.  Also  near  the  end  of  the  set,  HAVSTAR  1  was 
uploaded  with  the  result  that  the  last  several  distances 
were  very  large,  as  shown  in  figure  7.16. 

In  set  5,  a  very  high  ORE  value  for  NA7STAR  1 
contributed  to  the  extremely  large  position  difference 
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between  that  determined  by  the  theodolites  and  that  for  the 
differential  GPS  position  {fig.  7.47).  The  aziauth  remained 
nearly  constant  between  the  two  positions  with  a  median 
aziauth  of  43°  (fig.  7.48).  These  values  for  set  5  show 
that  the  large  ORE  for  HA7STAR  1  was  the  aost  important 
contributing  error  source  (fig.  7  .16)  .  The  large  tJRE  may 
have  been  caused  by  activation  of  the  hydrazine 
stabilization  jets  to  dampen  momentum  on  NAYSTAR  1. 

?or  set  6,  the  mean  distance  difference  between  the  two 
positions  was  60.0  aeters  and  a  median  aziauth  of  93°  (fig. 
7.49  and  7.50).  These  values  indicate  that  the  upload  of 
NAYSTAR  1  near  the  beginning  of  the  set  (0715  on  fig.  7.16) 
was  received  and  utilized  in  the  Hanpack  position  algorithm, 
but  the  higher  distance  between  the  positions  also  indicates 
a  large  UBE  (also  shown  on  fig.  7.16). 

Data  froa  set  7  were  not  evaluated  because  the 
stationary  Hanpack  receiver  did  not  record  the  GPS  position 
for  each  fix  due  to  radio  failure.  The  aean  distance 
between  the  two  positions  for  set  8  was  121  aeters  with  a 
aean  aziauth  of  68°  (fig.  7.51  and  7.52).  Both  the 
distances  and  azimuths  show  that  the  ORE  for  NAYSTAR  1  at 
0530  (fig.  7.17)  was  increasing  in  magnitude  and  thus 
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Figure  7.51 
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Figure  7.52 
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resulting  in  a  greater  positional  difference.  Another 
factor  vas  the  bias  introduced  by  the  Sanpack  24-second 
reception  period.  This  is  evident  in  the  azinuth  for  the 
theodolite  ship  position  obtained  from  the  GPS  shipboard 
receiver. 

The  data  for  sets  9  and  10  vere  combined  for  evaluation 
and  the  results  mere  a  mean  distance  between  positions  of 
241  meters  and  an  azimuth  of  267°  (fig.  7.53  and  7.54). 
Values  of  that  magnitude  show  that  the  UBE  for  NAVSTAB  1  vas 
the  error  source.  Data  vere  sampled  from  0600  to  0640  (fig. 
7.17).  Thus  the  correlation  between  the  high  HA7STAB  1  UBE 
and  the  high  mean  distance  difference  between  positions. 
Again,  receiver  bias  also  introduced  an  error  into  the 
evaluation. 

Analysis  of  set  11  shows  a  bimodal  distribution  in  both 
the  two  position  differences  and  azimuths  (fig*  7.55  and 
7.56) .  The  peaks  on  the  distance  histogram  are 
approximately  133  and  160  meters;  the  peaks  on  the  azinuth 
histogram,  approximately  220  0  and  236° .  Both  peaks  are 
caused  by  the  upload  of  NAVSTAB  3  halfway  through  the  set 
and  the  increasing  magnitude  of  the  NAVSTAB  1  BBS  just  after 
the  upload  as  shown  by  figure  7,17.  The  total  positional 
difference  for  set  11  was  approximately  154  meters. 
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Figure  7.53 
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Figure  7.56 
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Evaluation  of  the  combined  sets  12  and  13  (fig.  7.57  and 
7.58)  result  in  a  mean  distance  between  the  two  positions  of 
172  aeters  and  a  aean  aziauth  of  262° .  The  aean  distance 
was  large  because  of  the  large  aagnitude  of  the  NAVSTAR  1 
QBE  (fig.  7.17).  A  contributing  source  of  error  was  the 
bias  introduced  by  the  Sanpack  receiver. 

Overall,  the  data  using  a  five-satellite  constellation 
were  as  good  or  better  than  the  data  using  only  four 
satellites.  The  results  for  12  Say  using  five  satellites 
(fig.  7.59)  show  a  aean  distance  of  11.4  aeters  and  a  aedian 
of  10.6  aeters.  This  conpares  favorably  with  the  results 
from  sets  1  and  2  (fig.  7.39  and  7.40),  which  had  aean 
distances  of  10.1  and  11.7  aeters,  respectively.  The  five 
satellite  data  froa  13  Say  showed  a  very  large  variance  due 
to  the  large  NA7STAR  1  ORE.  Therefore,  data  were  not 
coapared  to  reaaining  13  Say  data.  Data  on  15  Say  for  the 
five  satellite  constellation  had  a  biaodal  distribution 
because  NAVSTAR  1  was  uploaded.  Once  the  uploading 
occurred,  the  peak  distance  was  137  aeters  (fig.  7.60)  which 
coapares  favorably  with  121,  241,  154,  and  172  aeters  for 
sets  8,  9  and  10,  11,  12  and  13.  Thus  with  a  five-satellite 
constellation,  positioning  results  were  as  good  or  better 
than  the  four- satellite  constellation. 
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Figure  7.57 
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Figure  7.58 
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Figure  7.59 
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Figures  7.61  to  7.74  show  the  difference  by  fix  between 
the  latitude  and  longitude  plots  of  the  theodolite  ship 
position,  the  differential  SP5  position,  and  the  uncorrected 
GPS  position.  During  aost  of  the  testing,  the  longitude  of 
both  eethods  of  GPS  positioning  was  in  closer  agreeaent  with 
the  theodolite  ship  position  than  was  the  latitude.  Shown 
on  figure  7.9  and  7.10  are  the  satellite  elevation  and 
aziauth  data  for  the  testing  period.  This  data  indicates 
that  the  satellite  constellation  had  a  larger  effect  on 
latitude  than  longitude  because  the  predoainant  ranging  was 
coaputed  froa  one  quadrant. 

The  large  differences  between  the  GPS  position  relative 
to  the  theodolite  ship  position  were  easily  identified 
during  the  shipboard  testing  period.  Latitude  and  longitude 
plots  (fig.  7.65,  7.66,  7.69  through  7.74)  for  sets  4,  5,  9, 
10,  and  12  show  the  deviations  of  the  GPS  receiver  positions 
froa  the  predicted  positions  that  were  obvious  to  the  GPS 
operators  during  the  shipboard  testing.  For  sets  4  and  5, 
the  stationary  GPS  receiver  positions  began  to  deviate 
greatly  froa  the  geoceiver  position,  resulting  in  a  large 
differential  correction  to  the  GPS  ship  position.  These 
large  differential  corrections  caused  the  differential  GPS 
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The  objective  of  the  study  vas  to  determine  whether  the 
application  of  GPS  in  a  translocated  or  differential  node 
provides  sufficient  position  accuracy  for  near-shore 
hydrographic  surveys.  Several  assumptions  were  aade  prior 
to  the  test:  (1)  two  GPS  receivers  required  for  the 
differential  aode  performed  the  computations  similarly,  so 
that  in  the  worst  case  a  constant  bias  would  be  displayed  if 
both  receivers  were  colocated;  (2)  the  process  of  updating 
the  epheaeris  information  would  not  contribute  a  significant 
positional  error  to  GPS  in  the  differential  aode;  (3) 
reasonable  accuracy  could  be  obtained  by  using  the  output 
froa  the  receiver  control  display  unit,  that  is, 
sophisticated  processing  net  hods  were  not  required;  and  (4) 
the  Nanpack  receiver  results  would  not  differ  when  used  in 
stationary  or  low-dynanic  environments  (shipboard  at  4 
knots) . 

None  of  these  assumptions  proved  to  be  valid.  It  was 
originally  expected  that  the  two  receivers  would  produce 
siailar  positions  if  both  sets  were  operating  at  the  saae 
tiae  over  the  saae  location.  One  limitation  of  these 
receivers,  however,  is  that  current  epheaeris  inforaation  is 
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not  received  unless  the  receiver  has  been  activated  prior  to 
the  update.  [Ref.  66]  For  a  period  of  tise  after  the 
receiver  is  activated,  the  portion  of  the  satellite  sessage 
that  contains  the  updated  epheseris  is  not  acknowledged. 
Therefore,  position  degrades  rapidly  when  a  flAVSTAR  is 
updated  and  the  current  epheseris  information  is  not 
utilized  by  the  receiver.  For  this  study,  unfortunately, 
the  initial  tiaes  of  activation  of  the  two  receivers  were 
not  recorded.  Hence,  the  effect  on  position  caused  by  this 
receiver  paraneter  was  not  defined. 

The  ORE  plots  illustrate  that  the  epheaeris  updates  did 
contribute  significantly  to  the  range  error  and, 
consequently,  the  resulting  position  error.  The  position 
determined  by  the  Manpack  receivers  was  not  affected  as 
seriously  by  the  clock  error  in  NAVSTAR  1  when  a 
five-satellite  constellation  was  used.  The  inaccuracy  of 
pseudo- ranges  received  froa  this  satellite  required  frequent 
updates  froa  Vandenburg,  especially  on  13  Hay.  A  receiver 
that  had  the  option  of  totally  eliainating  the  erroneous 
pseudo-range  would  have  given  the  best  results  for  the 
study.  Since  the  receivers  operated  independently,  and  the 
epheaeris  inforaation  was  not  concurrent,  the  differential 
aode  was  effective  for  only  a  few  cases. 
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k  siaple  correction  was  applied  to  the  values  output  on 
the  control  display  unit  of  the  Hanpack  receiver.  This  did 
not  per ait  application  of  tropospheric  or  refraction  error 
corrections,  other  than  that  applied  by  the  receiver  during 
coaputation  of  the  position.  k  sore  sophisticated 
processing  method,  resulting  in  greater  accuracy,  would 
require  acquisition  of  the  pseudo- ranges  and  subsequent 
postprocessing  to  deteraine  the  receiver  positions. 

The  best  results  during  the  onboard  portion  of  the  study 
were  obtained  when  the  ship  was  allowed  to  drift  at 
approxiaately  1  knot  on  12  Say. 

In  conclusion,  use  of  two  Phase  I  Hanpack  receivers  in 
the  differential  aode  aeets  the  required  accuracy  of  10 
aeter3  only  if  the  ship  speed  is  less  than  1  knot  and  other 
paraaeters  such  as  epheaeris  inforaation  and  the  nunber  of 
satellites  used  for  the  position  coaputation  introduce  a 
ainiaua  amount  of  error.  Further  testing  using  differential 
GPS  should  incorporate  aore  sophisticated  processing  aethods 
and  aore  advanced  receivers  to  achieve  the  survey  accuracy 
required  for  nearshore  hydrographic  surveys. 
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APPENDIX  A.  DESCRIPTION  OF  ERROR  THEORY 

Proa  application  of  error  theory  in  the  evaluation  of 
positional  inforaation,  it  is  possible  to  establish  a 
aeaningful  accuracy  state aent  subject  to  unifora 
interpretation.  To  provide  a  logical  and  acceptable  basis 
POB  COMPUTATION  AND  COMPARISON,  POSITIONAL  ERRORS  ARB 
assuaed  to  follow  a  noraal  distribution.  The  assuaption  is 
valid  because  positional  error  coaponents  generally  follow  a 
noraal  distribution  pattern  when  sufficient  data  is 
available. 

The  statistical  treataent  of  errors  utilizes  aeasurable 
errors  obtained  in  the  sources  of  positioning  inforaation. 
Analysis  of  the  linear  coaponents  provides  a  two-diaensional 
expression  of  the  accuracy  of  the  positioning  systea  in  the 
fora  of  an  error  ellipse.  The  use  of  an  error  ellipse  is 
coaplicated  by  the  problea  of  axes  orientation  and 
propagation  of  elliptical  errors.  Therefore,  the  ellipse  is 
coaaonly  replaced  by  circular  fora  which  is  easier  to  use 
and  understand.  The  linear  standard  errors  are  coabined  and 
converted  to  a  circular  distribution  and  the  final  accuracy 
stateaent  is  expressed  in  teras  of  circular  errors. 
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The  circular  fora  assuaes  the  circular  probability 


distribution  function.  This  function  expresses  the 

probability  that  the  radial  error  will  be  equal  to  or  less 
than  the  radius  of  the  probability  circle.  The  dispersion 
of  errors  within  the  distribution  is  aeasured  in  teras  of 
precision  indexes,  each  of  which  represents  the  error  which 
is  unlikely  to  be  exceeded  for  a  given  probability.  The 
preferred  circular  precision  indexes  are  the  circular 
standard  error  (a  )  ,  the  circular  probable  error  (CPE  or 
CEP),  the  circular  aap  accuracy  standard  (CNAS) ,  and  the 
circular  near- certainty  error,  or  3.5  sigma  [3.5  a  )  .  For 
each  precision  index,  the  given  probability  is  39.35%,  50%, 

90%,  and  99.78%,  respectively. 

Each  of  the  precision  indexes  can  be  converted  to  one  of 
the  others  by  using  the  following  table: 

Circular  Error  Conversion  Factors 


39.35% 

50% 

99.78% 

39.35% 

1.0000 

1.1774 

2.1460 

3.5000 

50% 

0.8493 

1.0000 

1.8227 

2.9726 

90% 

0.4660 

0.5486 

1.0000 

1.6309 

99.78%  0.2857  0.3364  0.6131  1.0000 


The  rapid  approximation  for  determining  the  circular 
standard  error  is 

CSE  =  0.5000  (ox  ♦  ay) 

where  a  x  and  v y  are  linear  standard  errors  obtained  froa  the 
sources  of  positioning  information,  and  therefore  are  the 
axes  of  the  error  ellipse.  For  a  truly  circular 
distribution,  where  0  c  is  identical  to  the  circular  standard 
error,  o  x  and o  y  are  equal  and  the  angle  of  intersection 
between  the  two  is  90°.  In  all  other  cases,  a  noraal 
circular  error  distribution  is  substituted  for  the 
elliptical  distribution.  The  substitution  is  satisfactory 
for  an  analysis  within  specifiedaain/oaax  ratios,  where o  ain 
is  the  minor  or  smaller  linear  standard  error  of  the  two. 
Distortion  in  the  error  distribution  between  elliptical  and 
circular  forms  at  lower  o  ain/  o  aax  ratios  Units  the  use  of 
the  circular  concept  to  ratios  greater  than  0.2  . 

The  CMAS  is  the  precision  index  used  as  the  U.S. 
National  Map  Accuracy  Standard.  It  is  interpreted  as 
limiting  the  size  of  error  which  90*  of  the  positions  will 
not  exceed.  CMAS  is  computed  from  circular  standard  error 
as  follows: 


CMAS  =  2. 1460  a  c 


The  two-diaensional  error  distributions  for  positions 
deterained  by  the  &RTE9IS  systea  and  the  theodolite  network 
are  described  by  CHAS.  The  accuracy  of  the  specific 
positions  is  expressed  by  a  stateaent  of  probability  and 
error  aagnitude.  The  accuracy  stateaent  does  not  aean  that 
the  error  in  position  is  exactly  the  value  shown  in  figure 
5.1,  rather  it  expresses  the  probability  that  the  error  in 
position  will  not  be  larger  than  the  error  given. 
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